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A SURVEY OF STABILITY ANALYSIS TECHNIQUES FOR 
AUTOMATICALLY CONTROLLED AIRCRAFT 
By Arthur L. Jones and Benjamin R. Briggs 


SUMMARY 


A survey of the stability analysis techniqxies for automatically con- 
trolled aircraft is presented. The survey is limited to the techniqvies 
commonly applied to linear, continuous-control systems wherein the dif- 
ference between the oul 5 )ut and input quantities is measured continuovisly 
and is used in the operation of the system (a closed— loop system) . An 
evalviation of the techniques, based on the kind and amoimt of information 
derivable, is included. An illustrative example is also presented to 
demonstrate the calculations involved for a typical aircraft— autopilot 
combination. 


INTRODUCTION 


The application of automatic control to the operation of aircraft 
has complicated the analysis of the stability of the aircraft motion. 
Before this innovation no atten^jt had been made to analyze the effect 
of directed control, as provided by the human pilot, on the aircraft 
stability. The omission of this important factor was due for the most 
part to the apparent impossibility of specifyi^ htmian response chsirac- 
teristics.^ Consequently, aircraft stability was evaliiated oh the basis 
of the \inregulated aircraft motion (or in the case of flight tests on 
the basis of pilots’ opinions). ' This type of analysis was carried out 
both with the controls fixed and with the controls free in order to 
provide as comprehensive a test of the stability as possible. 

The combination of the pilot and aircraft represents a closed loop 
of operations; that is, the pilot observes some output q\iantity of the 
aircraft such as its attitude, compsires it with the desired attitude, 
and operates the controls to reduce the difference or error between idle 
actual and desired attitudes. The airplane then changes attitude and 
the pilot repeats the process of observation, comparison, and control 
operation until the error is reduced to zero. If the human pilot is 
replaced by an ::utopilot, consisting essentially of an error-measuring 
^Attempts recently have been made to study man as an element of a con- 
trol system. (See references 1 and 2.) 
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device and a servo motor, the same closed loop of operations takes 
place. Furthermore, the response chsiracteristics of the autopilot can 
he measured or estimated and this information, together with the response 
characteristics of the aircraft, makes It possible for the autopilot— 
aircraft cosibination to be* analyzed as a closed— loop automatic— control 
system or, more concisely, as a servomechanism. In the literatiire 
relating to automatic controls there have been atten^jts made to distin- 
guish between servomechginisms and sin 5 >le automatic control systems often 
described as regulators. (See reference 3») The distinction is quite 
subtle, however, and has no particular significance in the analysis of 
automatically controlled aircraft. 

At the present time there are two methods of approach in the 
anal y sis of the combined aircraft-autopilot stability. One approach 
follows very closely the "equations of motion" type of analysis familieir 
to the aeronautical engineers and accomts for the autopilot by means of 
added stability derivatives, forcing fimctions, or supplementary eqvia— 
tions. (See references 4, and. 6 .) The other approach follows along 
the line developed by the mechanical and electrical engineers in their 
analyses of automatic- control systems and treats the aircraft as merely 
one component of the control loop. (See references 7, 8 , and 9 .) Both 
of these methods are considered herein, and an effort is made to cata- 
logue and correlate the various techniqvies that tave been developed. 

The techniques covered might well he described as elementary since in 
the actual process of design and construction of automatically controlled 
aircraft the analysis techniques often become more elaborate due to the 
extensions and modifications that are employed. 

This survey will be concerned with continuous-control systems only. 
The discontinuous or "onr-off" type of system may be designed to perform 
similar operations, but the analysis of such systems differs consider- 
ably from that of a continuous system. (See reference 10.) Moreover, 
only linear systems will be considered in this report. Nonlinearities 
complicate the analysis to the extent that one particular system woxild 
probably require as much time for analysis as a whole class of linear 
systems, * 

The report will be devoted primarily to discussions of the tech- 
niques of analysis of the stability and performance of automatically 
controlled aircraft in contrast to synthesis techniques for such systems. 
Analysis, being merely an evaluation of the existing characteristics of 
a system, is a simpler procedure than synthesis which is concerned with 
building up of components or the determination of modifications that 
will provide a system meeting certain standards of stability and per- 
formance. The definitions of such standards or criteria are a pressing 
problem in the development of automatically controlled aircraft at the 
present time. No attenpt will be made herein to formulate such standards 
but the fundamental parameters that probably WQVild be used in the for- 
mulations will be pointed out and the analyses will be judged on the 
basis of how completely and how readily these parameters may be evaluated 
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The general plan of the report vill he to present a brief descrip- 
tion of the conponents and analytical concepts associated vith servo- 
mechanisms followed hy a discussion of the stability analysis technlq,ues* 
i^ppendix A contains an application of the techniqiies to a typical 
aircraft-autopilot combination. The precise definitions of terms and 
the detailed descriptions of the techniques, however, will be left to 
the many eind varied references that are mentioned. The symbols 
coefficients used in the text are listed and defined ‘in appendix B. 

COMPONENTS AND CONCEPTS 

In the subsequent discussion the automatically controlled aircraft 
is considered to be a rather complicated servomechanism. There are many 
existing definitions of servomechanisms. The following definition suits 
the nomenclature enployed in this report and is based on those given in 
references 3# 8, and 11: A servomechanism is a combination of elements 

controlling a source of power which tends to make the output of the 
device follow a particular input by use of a command signal, the strength 
of ^ich is a function of the error between the input signal and the 
fed-back output signal. 

Many practical exasples of simple servomechanisms are discussed and 
illustrated in references 3, 11, 12, 13,14, and 15. For the purposes 
of this report the following block diagram, with all the couponents 
lumped into one box, will serve as the basic schematic • illustration of 
the servomechanisms discussed herein. 



6 symbol for the output variable 
d± symbol for the input variable 
E symbol for the error variable, ©i — 0 

The source of power mentioned in the preceding definition corre- 
sponds to the autopilot or to the servo motor part of the serVomecLunism. 
The sensing device for the production of an error signal is also part of 
the autopilot.- The airplane is the controlled -object and its motion is 
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the output. The control surfaces of the airplane, actuated by the auto- 
pilot, provide the aerodynamic control of the system. The input signal 
is usually derived from a fixed course that the aircraft is trying to 
follow or from a target that the aircraft is trying to intercept. 

The feeding back of the output to the sensing device provides a 
follow-^q) characteristic which makes the System a closed loop and thereby 
provides the form of control that is the basic characteristic of any 
servomechanism or regulator. In open— loop control systems, the control 
acts according to some preset command and is not at all governed by its 
output action. 

In general, servomechanisms may be designed to control quantities 
other than output position. For instance, the output velocity, a force, 
a press\rre, or almost any physical quantity can be the controlled 
variable. (See reference 3«) In application to aircraft) however, 
automatic control will refer almost exclusively to the position of the 
aircraft. Consequently, the aircraft-autopilot combination may be 
referred to as an error-sensitive, follow-tp, remote-control, position- 
ing servomechanism. 


PHYSICAL COMPOHENTS OF M AUTOMATICALLY COMEOLLED AIRCRAFT 


As previously indicated, all automatic control systems consist of 
a unit or "plant” to be controlled and a unit or "controller" to do the 
controlling. Tie boundaries separating these conponents are not too 
well defined for the many and varied mechanisms that could be classified 
as automatic control systems. For aircraft application, howevi.r*, it is 
probably most convenient to consider the aircraft to be the plant and 
the rest of the system to be the controller. The controller could be 
further separated into a measuring device and control imit consisting 
of an anplifier and a servo motor. A block diagrsun repjresenting the 
above breakdown of th.e automatically controlled aircraft is given as 
follows : 
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where 6 is the deflection of the aircraft control svirface as actxmted 
by the control imit. Not all of these components are essential to a 
servomechanism in general, nor are all servomechanisms restricted to 
these components, but the breakdown used will be suitable for the cases 
considered herein. 


Error-Measuring Devices 


The function of the erroMiieasiirlng device is to coD^iare the fed- 
back output signal with the input signal and thus determine the error 
between the desired position and the actual position of the aircraft, 

The controlled variable could be velocity rather than position but, 
for aircraft, control of the position is usually desired. 

Error-detecting devices can be mechanical, electrical, or of many 
other forms. (See references 3, 12, 13, l6, 17, l8, and I 9 .) In air- 
craft application, however, a gyroscopic device is most often used. The 
tendency of the axis of the rotor to remain fixed relative to an intial 
direction in space is utilized to indicate deviations from a set coiorse. 
Various other instruments such as anemometers, pendulums, etc., can be 
eind are en 5 >loyed. 


Control Unit 


The control unit in the autopilot-aircraft servomechanism transmits 
the error signal to the plant with suitable amplification and power to 
produce the desired output from the plant. This unit is in itself a 
servomechanism with its own inner feedback loop and usually consists of 
an eimplifier, lead or lag networks, and a servo motor. 

The servo motor, which can be. of electrical, hydraulic, or mec ha nical 
form, provides the torque to operate the controls of the airplane. The 
amplifier converts the relatively weak error signal to a CTirrent sviffi— 
cient to operate a hydraulic valve or electric device used to control 
the servo motor. 

Phase lead or lag networks may be incorporated within the inner 
loop to Improve the damping and steady-state error of the control-unit 
servomechanism or may also be added external to the inner loop to pro- 
vide definite characteristics of damping for the autopilot-airplane 
combination. The properties of such networks are 'discussed at greater 
length in later sections of this report. , 
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Controlled Plant 


Whether the aircraft he a conventional airplane, a powered missile 
capable of sustaining level flight, or a drop bomb with some means of 
guidance, the dynamic control is achieved by the deflection of aero- 
dynamic control stirfaces. Both flap— type control surfaces and variable- 
incidence control surfaces are \ised. With the exception of a variable- 
incidence lifting wing, these control surfaces are usxially mounted near 
the tips of the wings and near the extremities of the body in order to 
provide moments of sufficient magnitude. The control of the aircraft, 
therefore, is usually obtained by a torque producing one of three rota- 
tions: roll about the longitudinal axis using ailerons, pitch about the 

lateral axis using elevators, or yaw about the vertical axis using a 
i^idder or rudders. 

The aircraft with these control surfaces is a much more complicated 
plant to control than those usually considered in books and articles on 
servomechanisms. In particular, more than one of these controls can be 
in operation at the same time. Thus, tq> to three degrees of control, 
two lateral and one longitudinal, must be dealt with as well as possibly 
six degrees of freedom. An additional degree of control is involved if 
both an elevator and a variable-incidence wing are used in the longitu- 
dinal control. 

The combined effects of the actions of two controls, whether they 
are actuated by separate autopilots or by one autopilot, will complicate 
the analysis if the resiilting modes of motions are Interdependent. This 
situation can occ\ar in the lateral case, for instance, if the ailerons 
and rudder are simultaneously employed or in the case of a missile so 
rigged as to destroy the independence of lateral and longitudinal motions 
provided by the lateral symmetry of most aircraft. Analyses based on 
one degree of freedom and one degree of control can be extended in 
special cases, or if approximations axe employed to the case of auto- 
matically controlled aircraft, as will be evident from the subseqxient 
discvissions of the techniques Involved. Methods that have been advanced 
for coping directly with more than one degree of control, however, are 
rare and will be treated apart from the commonly considered one degree 
of control problem. 


FUHDAMEHTAL AHALTTICAL CONCEPTS 


The four dynamic response parameters that shall be considered as 
the basic items for the evaluation of stability and performance charac- 
teristics are (l) speed of response, (2) steady— state error, (3) time 
to damp, and (1+) amount of initial overshoot. Reference to these para- 
meters will help to clarify the functional purpose of the various 
components of the control loop in addition to providing a basis for the 
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future specification of stability and performance standards. These 
four items are best described in terms of the transient response (time 
history) of the motion, and they are illustrated qualitatively in the 
following diagram of a so-called "underdaraped" (references 3 and 13) 
oscillation: 



vhere 

^ approximately the speed of response 
‘'r 

Sx steady— state error 

tj ^^2 to damp to half amplitude 

Sg initial overshoot 

For aperiodic motion the overshoot is, of course, missing. The time to 
damp is often expressed in other ways such as the time for the amplit\ide 
to damp by l/e, the cycles to damp to l/2 or by l/e, etc. 

The first two dynamic response peirameters sure affected mainly by the 
type of control used in the servomechanism eind the last two by the type 
of damping employed, although none of the fovir parameters is completely 
independent of either control or damping. Also, both the stability and 
performance are related to all four parameters'. The stability of the 
system is, of course, almost directly related to the damping. It is 
impossible to have stability without somtie positive means of control, 
however, whether It be proportional, derivative, integral, or otherwise. 
The problem of stabilizing a servo system, consequently, is inseparably 
involved with considerations of both control and damping. The improve- 
ment of the perfor man ce of a servo system is sijnilarly involved. 

In the concept of stability used by the aerodynamic 1st,' any positive 
damping of all the modes of motion (totaling the same as the number of 
real and conjugate complex pairs of roots of the characteristic equation), 
oscillatory or aperiodic, constitutes stability. This corresponds to 
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having all negative real parts for the roots of the characteristic equa- 
tion of the system. This definition of stability is quite common. An 
exception to this definition was used by Whiteley in reference 20. 
Whiteley considers a system to he stable if the response to a distiirhance 
is nonoscillatory with small initial overshoot. The first definition 
seems to he the most sxiitahle and the most generally applied to stability 
and therefore will he employed throughout this report. Whiteley' s defi- 
nition is reaU.y a definition of satisfactory performance and as such 
co\Jld he usefTilly enployed. The establishment of suitable performance 
standards, however, is a controversial issue at the present time through- 
out servomechanism literat\ire, and the general feeling seems to he that 
the standards will have to fit the intent of the design. Again it shovild 
he einphasized that no attempt will he made in this report to establish 
standards for the stability and performance of automatically controlled 
aircraft. 

Time lags (usually of the exponential type in a continuous control 
mechanism) usually affect the stability of the system adversely, and, 
as one or more of these time lags are usually present in a servo- 
mechanism, the control and damping problem is immediately evident. In 
the following sections of the report, discussions of the various types 
of time lags and of the various types of control and danping will be 
presented. Inasmuch as pure derivative or integral control is impos- 
sible to obtain and the apparatus intended to provide such control is 
often complicated, discussions of passive lag or lead networks and of 
feedback networks that will give suitable approximations to integral or 
derivative control are included. 


Transfer Functions 


Althoxi^ stability analysis is concerned ultimately with the servo 
system as a whole, the investigation of the contribution of each conpp— 
nent in the loop to the stability is an important part of the analysis." 
In order to make such a study it is convenient to utilize the concept 
of the transfer function. 

A trsnsfer function can be thought of as a mathematical operator 
■vdiich relates an output to an input. The operator represents the action 
of a physical part of the servo loop and it can be obtained in operator 
form by taking the Laplace transform of the differential equation relat- 
ing the output to the input. Certain forms of the transfer function 
also can be obtained by analyzing experimental data. Mathematically 
speaking, there are many different foms of the transfer function in 
use in servomechanism analysis, and unfortunately there are some incon- 
sistencies in the terminology associated with these functions. Before 
discussing any of the differences in terminology, however, the two most 
common foms of the transfer function will be derived. 
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For perfectly arbitrary inputs, the differential eqmtion itself 
can be considered as an elementary type of transfer function. More pre- 
cisely, however, a transfer function is obtained by transforming the 
equation using either the Heaviside operator or the Laplace transform 
technique. The transfer function in this form is extremely usefiil 
because of its manageability. To illustrate the derivation of the oper- 
ator type of transfer function, consider a second— order servomechanism for 
which the following differential eqmtion holds 

0 (t) + 2 ^Uq 6(t) = E(t) ' ( 1 ) 

Assuming that the initial values are all zero, the Laplace transform of 
the eqviation is 

- (p2 + 2Cwn p) e(p) = (2) 

Now rewriting the transformed equation in the following form: 

w 2 

e(p)=-g E(p) (3) 

p + 2^UnP 

the function can be considered as operating on E(p) to pro— 

P + 2CwnP 

duce 0(p). This fiinction, therefore, is the transfer f\mction between 
E and 0. It is the open— loop transfer function because of the vari- 
ables involved. 

In this operator form the transfer function can be used to determine 
the effects of a fairly leirge and general class of inputs oh the. output 
of the system, (See reference 11.) The procedtrre is as follows: The 

input function is put in operator form, the necessary algebraic manlp — 
ulations are carried out, and the inverse of, the Laplace transform is 
taken to obtain the time^istory solution for the output. 

The second type of transfer fxmction oftentimes employed is the 
frequency-response form. It can be obtained directly from the operator 
form by substituting Iw for p. A formal mathematical derivation of 
this ftinction is given in reference 11, 

A few examples of the symbols used for the transfer ftmctions are 
given below in the frequency-response notation. These functions repre- 
sent, in particular, the opexjr-loqp response of a servomechanism but are 
typical for all applications of transfer functions: 

^ = KG(1w) (Reference 3) (^) 

^ = K |g(1u)| (Reference 3) 


(5) 
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E 

0 

E 

0 

E 

0 

E 

0 

E 

0 

E 

0 

E 

In any 
between the 


= B + iC 

(Reference 7) 

(6) 

= Re^^ 

(References 7 and 21) 

(7) 

= Y(iu) 

(Reference 11) 

( 8 ) 

= (fP§) 

(Reference 9) 

(9) 

= (AR)e^A 

(Reference 9) 

(10) 

= F(iw) 

(Reference 12) 

(11) 

= A(iu)) 


(12) 

case the symbols stand for an operator giving 
input and output of a system or a component of 

a relationship 
a system. 


In regard to the terminology associated with the transfer fiinctions, 
the soTorce making the greatest distinction between types is reference 11, 
Therein, the term "transfer function" is reserved for the Laplace trans- 
form form ^d the frequency— response form is called the "frequency- 
response function." In reference 3 the term "transfer function" is used 
for all applications except the output— over— input operator 9/91 for the 
closed loop. The latter is designated as the "frequency response" of 
the servomechanism. The terminology in reference 11 is similar to that 
in reference 3 with the exception that there is no partiCTilar designation 
other than "output-input response" given to the 9/9^ operator. In 
reference 12 the term "transference" is used for both the Laplace operator 
£ind frequency— response form of the transfer function. Reference 22 uses 
"system transfer function" for 9/9± and "loop transfer function" for 
0/E, 


The application of these functions to the analysis procedures will 
be discussed in detail in the section of this report on techniques of 
einalysis. 


Characteristic Equations 


The study of the characteristic eq\iation of a servomechanism is in 
itself a very convenient means of checking the stability of the system 
as will be shown in more detail in the discussion of the techniques of 
analysis. This equation can be formulated by writing the differential 
eq\iation describing the dynamic behavior of the complete closed— loop 
system, then reducing the homogeneous form of the equation to a 
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polynomial. In some cases such as when the lag operator to he descrihed 
subsequently Is used, the homogenous form reduces to a transcendental 
eqxiation. Generally the reduction is achieved either by substituting a 
solution of the form e^t into the eqimtion or by writing the equation in 
operator notation, assuming the initial conditions to be zero. The 
roots of the resulting algebraic or transcendental eqxiation can then be 
determined. (See references 3, 23, 24, and 25.) If the roots are real 
and negative, or if the real parts of the complex roots are negative, 
the system is known to be stable. 

The occurrence of a pair of conjugate complex roots in the solution 
of the characteristic eq\xation indicates an oscillatory mode of motion. 

A real root indicates an aperiodic mode. In addition to the damping 
information provided by the real part of the complex root, the frequency 
of the oscillatory mode can be determined from the imaginary part. 


Time Lags 


The discussion of time lags will be divided into two parts. The 
first part will be on delay time, the period during which no response 
to a signal takes place. The second part will cover time constants 
which relate to the lags caused by exponential response to a signal. 


Delay Time.- The study of the dynamics of mechanical systems is 
usually based on a differential equation that ideally represents the 
application of forces and the resultant acceleration of the masses or 
inertias. In some cases where the system involves a chain of events, 
however, an accurate estimate of the response of the system to a signal 
is not analytically representable unless certain delay times be con- 
sidered that are not taken into account in the basic differential equa- 
tion. (See references 26, 27, 28, and 29.) Delay time can be defined 
as the period between the start of a signal and the start of the result- 
ant output. In electrical systems this kind of time lag is usually 
enough to be insignificant, but in purely mechanical, hydraulic, or 
pneumatic systems it is occasionally necessary to evaluate the effects 
of delay times or dead times as they are sometimes called. 

There are a nxmiber of ways to incorporate delay time in the analysis 
of a servo system. The best way, of course, would he to write a dif- 
ferential equation or a set of simxiltaneous differential equations which 
would represent perfectly the characteristics of the system. Even a good 
approximation to this ideal situation, however, is possible only for the 
case of electrical or electrical-mechanical systems. The combined air- 
craft and autopilot system is one that usually can be treated as an 
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electrical-ffleclaanical system. If a frequency-response type of analysis 
is being used and the experimental results are available for both the 
aircraft and autopilot, these results can be combined to yield almost 
exactly the characteristics of the behavior of the coxapled conq)onents. 
Furthermore, it is fairly common practice to take the linear theoreti- 
cal expressions for the aircraft response and the autopilot response, 
combine them, and proceed with the analysis on the assumption that the 
time lags are either closely approximated by these expressions or suffi- 
ciently small to insure reasonably accurate results. 

In cases where an expression for some component of the system or 
some operation in the system cannot be written, even approximately, a 
finite time delay is often assumed and a term expressing this delay 
inserted into the differential eqxiation for the system (making it a 
differential-difference equation) . To illustrate this operation, con- 
sider a system defined by a differential eqviation of the form 

9 (t) = f(D) E(t) (13) 

If the initial conditions are assumed to be zero, the Laplace transform 
of this equation is 


L[q<t)] =.A(p) L[E(t)] ( 14 ) 

Now if there is a delay of seconds in the application of the 

operations represented by f(D), E(t — t,j^) must be substituted for E(t) 
in the above equations. Then by the translation theorem of Laplace 
transforms 

L [E(t - T^)] = L[E(t)] , ( 15 ) 

and equation (l 4 ) can be rewritten as 

L[cp(t)] = A(p) L[E(t -T<i) ] 

= A(p) e-^'^a L [E(t) ] (16) 

to take into account the time delay. The operator is often 

referred to as the lag operator. Unfortunately, the application of this 
operator yields a transcendental characteristic eqiiation and-, in general, 
complicates the analysis to such an extent that many investigators have 
tried using a three— tenn approximation to the power series expansion of 
this operator. (See references 6 and 30 .) Such an approximation could 
lead and, in some cases, was reported to lead to erroneous results. 

Proper applications of the lag operator are demonstrated in refer- 
ences 23 and 27. 
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Time Constants .— In contrast to the delay time, which represents a 
fixed time delay between the start of a signal and the resulting 
response, the time constant concept arises from a certain exponential 
time response to an input signal, that is, a sort of sluggishness of 
response that is characteristic of systems having inertia, inductance, 
or the like. More specifically, the speed of response of a member in 
the system often cein be related to a certain constant or constants that 
appear as coefficients in the transfer function of the member. These 
constants, for dimensional reasons, are known as time constants. 


For an illustration of a simple case, consider a member of the 
control loop that has a transfer function of first order 

' riiT • ■ 

The symbol T represents the time constant and,^ in general, it can be 
related to some physical parameter of the control— loop member, (See 
references 3 and 12,) The solution of the differential equations for 
this loop would contain this constant as the inverse coefficient of an 
exponential term in the manner 

A(t) = (18) , 

Thxis the magnitude of the time constant is seen to give an indication 
of either the sluggishness of the response or the rate of damping of 
the response, depending on which phase the motion is in — the approach 
to the desired condition or the settling down of the overshoot. 

It can be seen from equation (l8) that the time constant is related 
to the damping constant associated with oscillations arising in a second- 
order or higher system. To clarify this relationship consider a simple 
second-order closed— loop system defined by the eqiiation 

0*(t) + 2t(u^0(t) + u^20(t) =u^ae^(t) (19) 

The term 0(t) is the damping term and 25un is called the danping 
coefficient for reasons more fully explained in the section of this 
report on danping. For a second— order mechanical system it can be 
demonstrated that 2^^ is dimensionally eqvial to l/sec and that the 
solution to the equation of motion (I 9 ) is of the form 


0 = e cos (u^t + e ) 


( 20 ) 
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Comparison of the exponential coefficients of equations (I8) and (20) 
shows that the time constant T corresponds^ in this case, to twice the 
reciprocal of the damping coefficient. In many instances where higher- 
order systems are considered, it can be shown directly that the damping 
coefficients are functions of the time constants of the first—, second—, 
or third— order terms in the transfer functions. As the terms become of 
high order, however, the mathematical and Inferential relationships 
between the time constants of the transfer functions and the damping 
coefficients and other parameters characterizing the behavior of the 
complete system become more and more indistinct and the value of the 
"time constant" concept is lessened. 

Various interpretations of the significance of the time constant 
can be fotind in references 3, 12, aud 26. In reference 12 much attention 
is devoted to the effects of both sluggishness (corresponding to time 
constants) and' dead time (corresponding to delay time). The combined 
effects of these two phenomena are studied as well as their Individual 
contributions . 


Control Gearing 


The concept of control gearing is used primarily by the aeronauti- 
cal engineer. (See references 6, 21. j It is the ratio of static 

control— surface deflection, caused by the autopilot in response to an 
error signal, to the static error signal. It can be. expressed as 



5(i(o) 

E(io)) 


03 =0 


( 21 ) 


This parameter merely expresses the ratio of magnitudes of the input and 
output of the servo controller (autopilot) and. does not give any indi- 
cation of p'hase lag or lead. In terms of the transfer fvtnction of the 
autopilot it is 


k = 




( 22 ) 


The control gearing concept is especially useful in defining 
stability boundaries as was done in reference 6, which typifies the aero- 
nautical dynamic— stability approach to the problem of analyzing automati- 
cally controlled aircraft. It also has been employed in the frequency- 
response type of analysis in reference 7 where the critical control 
gearing kcr> which corresponds to neutral stability, was investigated. 


It should be pointed out that the sensitivity or gain, defined by . 
the symbol K in reference is not necessarily eqxial to the control 
gearing but is related to it by 
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K = 


I^L=o 


(23) 


vhere |g | is the absolute magnitude of the frequency dependent part of 
the transfer fionction of the autopilot es^ressed as KG(icn). 


In reference 31 > the inverse of the control gearing is called the 
static follow-^q) ratio and is used as a parameter. 


Types of Control 


In servomechanism tenninology the term "control” refers to the type 
of signal fed to the control unit. The intent of the signal, of course, 
is to produce an output that will coincide with the input. Consequently, 
the signal is generaJJ.y made a function of the error between the input 
and output , although control signals based on the input signal or dir— ■ 
ectly on the output signal have been -investigated and used. 

Ihe four kinds of control that are most commonly considered in 
analysis are: (l) proportional control; (2) proportional plus integral 

control; (3) proportional plus derivative control; and (4) proportional 
plus a combination of integral, and derivative control. (See refer- 
ences 13 and 20.) Discussions of more complicated forms of control can 
be found in references 20 and 32 . 

The use of proportional control in 5 )lies that a signal directly 
proportional to the error is controlling the source of power so that the 
control surface of the aircraft will be deflected in a direction tending 
to reduce the error. This results in zero steady— state error for a step 
input in the absence of external stiffness load. Addition of integral 
control, on the other hand, provides a zero steady— state error for a 
constant velocity input and also for a step input in the presence of 
external stiffness load. Integral control cannot be added to an unstable 
system as a direct means of stabilization, but indirectly it may be used 
to increase the stability of an already stable system by permitting the 
use o:^ a lower gain factor or control gearing. Addition of derivative 
control or error-rate control, as it is sometimes called, provides antic- 
ipation of the change in error that is to take place by feeding the 
derivative of the error signal to the power sotirce. This speeding up of 
the response by use of derivative control, in contrast to the slowing 
down caused by the use of integral control, improves the damping of the 
transient oscillations. The effect of adding both the derivative and 
integral control signals can best be expressed in the words of Whiteley 
(reference 20), "...it may be deduced that a combination integral and 
differential— of— error terms is a desirable arrangement, the former to 
control rates of response and achieve specified steady— state errors, and 
the latter to promote stability." 



l6 


MCA TU 2275 


The effects of these various types of control on the danping will 
he discussed in the following section on damping. 

The more subtle characteristics of these controls are discussed at 
some length in references 3^ 12, 13, 20, and 26. Methods of approxi- 

mating these controls by use of l6ig or lead networks or feedback loops 
are discussed in later sections of this report. 


Danping 


The danping of a servo system is very closely related to the 
stability of the system. Without danping of some kind, the oscillations 
usually arising as one of the modes of motion wovild be sustained if not 
divergent. The aperiodic modes of motion are also affected by the 
damping and care must be taken not to have too much danping present or 
these modes and the over— all response of the servo system may become 
sluggish. 

The terminology associated with damping seems to have been developed 
in the study of staple linear, second— order systems having viscous daop— 
ing. Some of the constants involved have been labeled with an assortment 
of names by various authors. In order to niaiie the following discussions 
on the damping of ordinary second-order servomechanisms and higher- 
ordered aircraft-autopilot systems as clear as possible, a short summary 
of the terminology used will be given. 


Terminology .— The differential equation for a linear second— order 
system may be written in the following form: 


d% „ dx ^ /r^l \ 

d — 5 + f ^ + |ix = 0 (24) 

dt 

where 

X some dependent variable 
t time, the independent variable 
J mass or moment of inertia 
f the damping coefficient 


the spring constant 
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The roots of the characteristic equation of the preceding differen- 
tial equation are 


-f ± 






2 - 


2 J 


(25) 


Two parameters can be chosen such that these roots can be expressed in 
terms of two quantities instead of three. These two parameters are u)n 
the natviral frequency of the system, and ^ a dimensionless damping 
ratio. In terms of the previous constants these parameters are 



(26) 


f _ f damping coefficient 

^cr ~ ^ critical damping coefficient 


(27) 


The critical damping coefficient fg^ is the value of f for which the 
two roots rj^ and rg become real and equal and the oscillations are 
changed into aperiodic modes. 


The differential equation can now be written in the following form: 


d^x 
dt ^ 


+ — + (OjjSx _ 


0 


(28) 


and the roots of the corresponding characteristic eqioation become 


= -^<Dn±n% y 1-^^ (29) 


where con^l— is the actual frequency oi of the damped oscillations 
or, as it is sometimes called, the freq\iency of the transient response. 
When ^ is small cOn is a good approximation to the actioal frequency. 
The equation for the time response is of the form 



x(t) = e^i't + Rg e^2t 

(30) 

or 

x(t) = e“ ^"^^(Rgcos 0Jn7l-^2t + 



R4 sin cOnyi-^2 

(31) 


The term e~ contributes the damping and the coefficient ^tDjj is 

called diversely the damping constant, the damping coefficient (which 
unfortunately is also the name associated with f in the differential 
equation) in reference 33 , the damping rate in reference 13, and occa- 
sionally the damping factor. 
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The inverse of this quantity I/^uJq can be referred to as a time 
constant. It represents the time reqviired for the amplitude of- the 
oscillations of the dependent variable to decrease by a factor of l/e. 

In this sense, it is similar to the time constants of the 'first-order 
terms discussed in the section on time constants. The number of cycles 
for t he am plitude to dan^) by a factor of l/e would be given by u/2n5ojjj 
or vl— 5^/2it^. The aeronautical engineer often uses the time to damp 
to 1/2 an^jlitude tj ^/2 as a damping criterion. This quantity is 
related to' in the following manner; 


_ 0.693 

1/2 ^(Dn 


(32) 


By taking the quotient of the period P -vAiere 



(33) 


and the time constant l/Ccon it is possible to obtain a dimensionless 
measure' of the dan5)ing called the logarithmic decrement 




2n^03n _ 

o> f/l=F 


It can be shown that 


d = log 


x(t) 

x(t+P) 


(3M 

(35) 


Hence, the logarithmic decrement is the logarithm of the ratio of dis- 
placements one period apart. The inverse of the logarithmic decrement 
l/d is equal to the cycles required for the amplitude to decrease to 
l/e of its value. 


Types of Dampin g.— There are two kinds of damping usually associ- 
ated with simple second— order servomechanisms, namely, viscous (iamp ing 
which is always present to some degree whether desired or not arui the 
so-called error-rate danq)ing. The latter is also known as derivative 
and kinetic damping. In some instances the viscous (^aTtipibg performs its 
stabilizing service without dissipating as heat loses too much of the 
power provided by the servo motor. When this favorable condition does 
not exist, however, there are two possible ways to reduce the energy 
losses of the viscous damping. One way is to use simulated viscous 
damping often refeired to as tachometric feedback, which absorbs no 
energy from the servo motor, to replace the viscous danping to as great 
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an extent as possible. For a detailed discussion of simulated viscous 
damping see reference I3. The other way to reduce the losses is to keep 
the viscous friction as small as possible and yet stabilize the system 
by supplementing the viscous friction with error-rate damping. 


The types of dan5)ing that are present in an aircraft-autopilot 
combination are shown in the following diagram: 

Total Damping 
of the Combination 


Autopilot 

Damping 


Viscous Error 

I Rate 


Mechanical Simulated 

The autopilot dan5)ing is similar to the damping of a simple second-order 
servomechanism. The aerodynamic damping corresponds closely to the ' 
actual viscous danping in the autopilot. When the aircraft and auto- 
pilot are combined to form a servomechanism, the resulting modes of 
motion and the damping of these inodes are. not simply the sum of those 
of the individual components. The feedback arrangement that is the basic 
element in obtaining automatic control changes the characteristics of the 
modes of motion although the total number of modes of motion is not 
changed. 

Some calculations showing the relative 'magnitudes of the dnTnping 
factors involved will be imdertaken to indicate these effects more 
explicitly. It will be advantageous to simplify the expression for the 
airplane — autopilot response in order to obtain a form that can be com- 
pared more or less directly to the second-order form of servomechanism 
previously discussed. The first step in this development is to comibine 
the three eq.uations of longitudinal motion. used in the illustrative 
example in appendix A with the assumed autopilot equation 

(0.058 D® + D + 37.5) 6 = -37.5 kE (36) 

■idiere 

k control gearing 
5 control surface deflection 



Aircraft 

Damping 

I 

Aerodynamic 
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This is easily done by inserting the autopilot equation as a forcing 
function in the equation for the angular degree of freedom (a stability 
analysis approach that will be described in the discussion of the tech- 
niques of stability analysis). The resulting forma of the equation are 


cos Oo - - 2 tDo^ 2 tD 0 + (^-20^^-205,^ sin - |^)u' 

= O(Lift) ( 37 ) 




) ^ (^Lo 

T y 


+ Cip^ sin ocq ^ ^ ^ 


T D|i ' =0 


(Drag) (38) 


Be 


m 


Ba 


a + 


BCjn 8Ctji 1 BC]H ^ 00 

Sa he he J 


+ u' 

Bu* 


^Cm / -j^ 37 -^ E 


S 6 V 0.058 + D + 37.5 


(Pitching Moment) ( 39 ) 


Now it is a very common procediire in the case of longitudinal motion in 
level flight to assume that the longitudinal perturbation velocity u' 
can be neglected and that conditions can be further controlled such that 
eqioation (38) can be eliminated. This procedure eliminates one degree 
of freedom and leaves the two following simultaneous differential equa- 
tions : • 


(- 


-Cip cos Oo 



a — 2t Da + 2t D0 = 0 


(HO) 


BCm 

a + 

Ba 


Ba 


m 


Da + 


Ba 



BCni 

Be ’ 



t^D^ 


k E 

S 6 \ 0.058 D^ + D + 37.5 


(Hi) 
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Now combining these to get one equation in 9 and E yields 

I? - «•) •’••■1.- S' • (% »• (i 5? - -)• 


2 aCm 


+ T D 


aCm / bCL\ 1 SCm 

2 + ( cos Oo + ) r- 


]} 


( 

V ^ 




^ ^Cm / k 37.5 E \ 

—{ — ] (42) 

/ • ^6 \ 0.058 D + D + 37.5/ 


In order to put this in the form of a torque balance equation as commonly 
written for the second— order servomechanism tD can be factored out of 
the left-hand side of the equation and the sighs of al 1 the terms reversed 
(to make the numerical values of the coefficients positive) : 

^ I (- 1) ^ cos ^ 


K - S=) J ^ - S') -s oo .g)(- i gs) . 


(cT„.oosa„-^^^2TD) (- 


Sc 


m 


/ 


k E 

\ 0.058 


37.5 ^ 

+ D + 37.5 / 


m 


In this form the coefficient of the angular velocity D0 is easily 
recognized as the aircraft equivalent of the damping coefficient of the 
second-order servomechanism equation. The autopilot damping is shown 
by the characteristic eqiiation of the autopilot in the denominator of 
the right-hand side of equation (43). In order to get some idea of the 
relative amounts of the danping that the airplane and autopilot possess 
uncombined, the numerical values used in the illustrative example of 
appendix A are inserted in equation (43) to give 
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[(0.0231) +(0.0342 +0.0656 +0.0704) (0.967) D + (1.278 + 0.187)] 9 

_ (5.325 + 1.933 D) (0.6299) 37.^ k E 
(0.967) d( 0.058 d^ + d + 37.5) . 


or 


[(0.0231) D^+(0.l645) D + 1.465] 9 = 


(5.325 + 1.933 D) (0.6299) 37.^ kE 
(0.967) ,D(0.058 + d + 37.5) 


Putting the eq.mtion in the desired form, 


(45) 


[(D^ + 7.12 D + 63.4)] 0 = 3T*? . 1 ‘ 933 , p) (0.6299) k E 

(0.0231) D (D^ + 17.25 D + 647) (0.058) (0.967) 

(46) 


and then in transfer function notation yields 


0 = AOg Ap E 

•where for the aiiplane' 

. ^ = ' (^.325 + 1.933 D) (0.6299) 

6 (0.967) D (D^ + 7.12 D + 63.4) (0.0231) 

yielding 

^‘“n = ^ = 3.56; = 8.0 


and where for the au'topilot 


yielding 


A_ = £ 37.^ k 

• E (0.058)(d®+ 17.25 D + 647) 


(47) 

(48) 


(49) 


= = 8.63; “n = = 25.5 

This coniparlson shows that the airplane has about half the damping con- 
stant that the autopilot has and has a natural frequency of 8.0 radians 
per second; whereas the autopilot has a nat\iral frequency of 25.5 radians 
per second, (The actual frequencies of the oscillatory modes are 7.2 and 
23.9 radians, respectively.) In terms of the time to damp to half ampli- 
tude, the airplane requires 0.I9 second; whereas the autopilot reqiiires 
only 0.08 second. 
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When the two units are combined into a complete automatic control 
system the behavior of the resulting servomechanism is not the sinrple 
svim of the behaviors of the individ\ial parts as- mentioned previously. 
The basic feedback loop to the error detector makes the equation for 
the response (references 3, 22, eind 34) 


e _ 9/E _ -^s -^p 
0i 1+0/E 1 + Ap 

The new characteristic equation is 1 + Ap = 0 which yields a 
fifth-order equation that is no longer factorable into distinct air- 
plane and autopilot functions. Thus the new modes of motion are a 
conplex combination of the autopilot and airplane modes. 


. Niamerically the new characteristic equation is 

(0,967) D [(0.0231) +(o.i 645 ) d + 1.465 ] (0.058 + d + 37.5) + 

(5.325 + 1.933 D) (0.6299) k 37.5 = 0 (51) 


or 


d 5 + 24.4 + 833 D® + 5700 + 76,200 kD + 97,100 k = 0 (52) 

The roots to this equation were determined for k = 1.0 and k = 0.2: 
k=l.o, ki=-l.39; ^2,3=-1*97 ± 10.5 ij ^4^5=-9.52 ± 22.8 i 
k=0.2, 4.28^ X.g^g=— 1.76 ± 1.89 k^^g=~8,29 i 24,8 i 

It can be seen from an inspection of the values for the dandling constants 
(the negative of the real parts of the roots) and the act ual frequencies 
(the numerical values of the imaginary parts of the roots) that the 
oscillatory mode corresponding to 14^5 has a close correspondence to 
the mode of oscillation of the autopilot alone. It is also the mode 
least affected by the variation in the strength of the feedback of the 
closed loop as shown by small changes in the danping and frequency char- 
acteristics caused by the variation in the control gearing k. The 
other oscillatory mode X.2^3 corresponds somewhat to the airplane mode 
although it is apparent that these two modes show large differences in 
the damping and frequency characteristics. The 3 mode and the 
aperiodic mode are quite markedly affected by the variation in 

feedback. It is important to note that with regard to the determination 
of the mode or modes that will be predominant in the output motion of 
the aircraft— autopilot combination, the magnitudes of the various modes 
must be calculated in addition to the damping and frequency characteris- 
tics. 
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A check on these results for k=l was made hy calculating the roots 
of the characteristic eq,uations resulting from the combination of the 
autopilot with the unmodified equations of motion for the aircraft and 
with the equations of motion modified by an approximate factorization of 
the stability quart ic. The values of the damping constants obtained eire ' 
given in table I. These additional analyses both revealed that in the 
uncombined condition the aircraft possessed a slowly divergent aperiodic 
mode. This divergent mode disappeared, however, when the autopilot arn ^ 
aircraft were combined in closed— loop operation. The values of the 
damping constants for the oscillatory modes were in good agreement among 
the three, analyses. 


Passive Networks in Cascade 


In cases where it is not suitable mechanically to add control in 
the form of derivatives or integrals of the error, the damping or per- 
formance of a system can be improved by the insertion of passive networks 
in cascade. Passive networks in cascade are networks containing no source 
or sink of energy that are inserted in series with the other elements of 
the loop. These networks are usmlly inserted near the input end of the 
main servo sequence and are referred to as "input networks" in refer- 
ence 20. These passive networks only approximate the effects of pure 
derivative— error or integral-error control by the addition of phase lead 
or lag to the system. The fact that such networks can be employed 
advantageously is emphasized in references 3 nnd 20. 

The process of analyzing the changes in a servomechanism loop 
required to obtain the proper system behavior can be done quite simply 
on the polar type of plots of the frequency— response methods of analysis 
(described under the section Techniques of Analysis). The transfer locus 
of the servomechanism must conform to a rather specific pattern in order 
that the system be stable and othejrwlse satisfactory. The shape of the 
curves can be changed over certain frequency ranges by the use of these 
passive networks. Thus a passive network is usually designed to give 
phase lag or lead in such a manner that when its transfer locus is added 
to that of the rest of the servo loop the resulting locus has the desired 
shape. In reference 3, this procedure is referred to as reshaping the G 
locus. To illustrate, consider a simple servomechanism with a frequency- 
response function of the following form: 


2. = A = ^ 

E ^ io)(iaffi+l) 


( 5 *^) 
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where is a time constant. The polar 

plot of this transfer fxznction is given in 
the sketch. As it stands, the locus indi- 
cates a zero position* error type of servo- 
mechanism because at the low frequencies 
it is asymptotic to the negative imaginary 
axis . Now if it were desired to change 
this servo system to a zero-velocity— error 
type> for which the low-frequency asymptote 
lies along the negative real axis, the fol- 
lowing transfer function corresponding to 
an integrating network could be used: 


A 2 = 1 + 

io) Tg 


(55) 


where Tg is a time constant of the net- 
work. The network locus has the shape 
shown in the accompanying sketch. Adding 
the network 'to the initial servo loop gives 
the following transfer fvmctionj 



- = A3A2 = ^ 

E T2(1 caf (idTi + 1) 


(56) 


and the transfer locus takes the shape 
indicated in the accompanying sketch. 


Feedback Networks 



When passive networks in cascade cannot be used to provide the 
approximate derivative or integral control desired, feedback networks 
are often used. (See reference 3, pp. 206, 2I7, 224. See also refer- 
ence 20.) In fact, it is frequently more desirable in practice to use 
feedback rather than cascaded networks. (See references 3, 20, 32.) 

The servomechanism Itself is a degenerative feedback system in^that the 
output in some form is fed back to reduce the error. An auxiliary 
feedback loop in the system can be either degenerative or regenerative 
although the latter requires a more critical adjustment. Feedback is 
usually obtained from a parallel circuit in the loop that provides for 
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the addition of a later operative q.'uantity to an eeurlier one, the earlier 
quantity usually being the error signal. Feedback in a broad sense can 
also be obtained by inserting a dynamic element in the basic servo- 
mechanism feedback from the output to the error-measuring device. (See 
reference 3, pp. 22k and 225.) Tb© parallel feedback circuits can be 
passive or dynamic although passive networks are usually preferred for 
the sake of their simplicity. 

Designing feedback networks to perform a particular control or 
stabilizing fimction or to modify the performance of the servo system 
in some particular way is not as straightforward and simple a procedure 
as was the derivation of passive lag or lead networks. Such a synthesis 
problem is most easily solved if the feedback network can be made sub- 
stantially equivalent to a suitable cascaded network. As an example of 
this process, consider a servo system of the following form: 




(58) 


If a p^allel feedback loop is added in the following manner 
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the openr-loop and closed— loop responses become 


9 

E 


l+AjAg 


( 59 ) 


and 


9 _ Ai 
9^ 1+A^Aq+ Ai 


(60) 


Now it is desired to find the form of A^ in terms of the equivalent 
modification to the basic loop by a cascaded network A„ shown as 
follows ; 



For the cascaded network the open^ and closed— loop responses are 


9 

E 


A1A3 


and 


_ A3A3 
l+A^Ag 


(61) 


(62) 


Thus to get As in terms of A3 the open-loop responses can be equated 
yielding 


A A 


Ai 

l+A 2^2 


(63) 


or 
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A 


2 


Ai-A^Ag 

A ^ A 

Al Ag 


( 64 ) 


Ag 




(65) 


for the feedback form as a function of the equivalent passive network. 

Whiteley in reference 20 says "Comhinations of input networks and 
feedback networks can of course be used; practical considerations may 
well determine what role shall be assigned to the input and feedback 
meshes. Generally it would be preferable to assign the differentiating 
functions to the input network, and the integrating functions to the 
feedback." 

A table of equivalent feedback and passive networks is given in 
reference 20. 


T H E 0 E Y 
TECmriQUES OF AMLYSIS 


In describing the techniques of stability analysis that are commonly 
applied to servomechanisms, it is convenient to establish at the outset 
the number of degrees of freedom and degrees of control to be investigated 
and to classify the techniques accordingly. The number of degrees of 
freedom of a system is equal to the number of independent variables 
required to describe the behavior of the system. The nimiber of degrees 
of control, on the other hand, is equal to the number of controlling 
elements that directly produce the output of the system. For instance, 
the conventional aircraft has three degrees of control provided by the 
rudder, elevator, and aileron. 

The analysis of ordinary second-order seryomechanisms, as generally 
carried out in the existing literature on the subject, is usually con- 
cerned with one degree of freedom and one degree of control; whereas the 
conventional aircraft has six degrees of freedom and three degrees of 
control. Fortunately, however, the motion of the conventional aircraft 
can be separated into two categories: longitudinal motion and lateral 
motion, by virtue of the lateral symmetry of the aircraft. (See refer- 
ence 35« ) The longitudinal motion generally is considered to have three 
degrees of freedom, a longitudinal velocity, a vertical velocity, and 
an angular velocity of pitch about the' lateral axis. In some instances' 
it is convenient and permissible to assume the longitudinal velocity 
constant and consider only two degrees of freedom. In the longitudinal 
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case the elevator usually provides the only degree of control. The 
lateral motion has three degrees of freedom, a lateral or sideslipping 
velocity, an angular velocity of roll about the longitudinal axis, and 
an angular velocity of yaw about the vertical axis. There are usually . 
two de^ees of control available in the lateral case: An aileron pro- 

duces rolling motion, and a rudder produces yawing motion. Analyses of 
the lateral motion of aircraft have been carried out based on the assump- 
tion that the rolling of yawing degree of freedom can be isolated, but in 
• general all the lateral motions are coupled to the extent that separation . 
is not advisable. 

As far as the classification of stability analysis techniques for 
automatically controlled aircraft is concerned, the number of degrees 
of freedom will not be used directly inasmuch as the linearized 
differential— equation type of analysis treats the aircraft as a whole so 
that all three degrees of freedom in one category are stabilized simul- 
taneously. Thus the discussion of the techniques of analysis will be 
divided into two parts: the first relating to one degree of control and 
the second to two degrees of control. As yet an aircraft with more than 
two degrees of control for one category of motion has not been produced. 

In the development of guided missiles, many unconventional designs 
have been evolved, some of which have more and some less geometrical 
symmetry than the lateral symmetry characteristic of conventional air- 
craft. In such cases the degrees of freedom, their coupling, and their 
relationship to the degrees of control must be established before the • 
analysis is attempted. It seems likely, however, that the classification 
of analysis techniques for such aircraft should be made on the basis of 
the numbers of degrees of control for any one category of motion. 

The autopilot is not always introduced in the form of a differential 
equation based on its physical and dynamic characteristics. When satis- 
factory mathematical expressions for the autopilot are not available, 
many investigators make use of a lag factor or a lag operator to simulate 
the action of the autopilot. These artifices were discussed to some 
extent in the section on delay time. It should be pointed out, however, 
that the use of the lag operator makes the characteristic equation , 
transcendental and consequently difficult to handle. Most of the tech- 
niques to be described can, in principle, handle the transcendental 
characteristic equation though it may not be practical to do so. 

In appendix A sample calculations are presented for a high— subsonic- 
speed fighter type of aircraft that include an application of each of 
the following techniques. 
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One Degree of Control 


In the consideration of automatic control it was only natural that 
the aeronautical engineer should attempt to perform an analysis by 
modifying the well-known set of simultaneous differential equations 
representing the motion of an airplane. (See references 35 and 36.) 

Two artifices were used. One was to include the autopilot action in the 
form of an added stability derivative (references 4 and 5 ) and the other 
was to include the autopilot action as part of a forcing fmction. (See 
reference 6.) These operations are equivalent. The obvious effect of 
adding the autopilot in this manner is to raise the order of the equa- 
tions and thereby introduce new modes of motion of the aircraft. 

A different approach is used by the servomechanism engineers, one 
that follows more or less the techniques of network or circuit analysis 
used to a great extent by the electrical engineers. In this approach 
the airplane is included as one member of the control loop and its dyna- 
mic response characteristics are lumped with those of the other compo- 
nents of the loop when the over-all response characteristics of the 
system are to be evaluated. 

As should be expected these two approaches provide essentially the • 
same initial equitions representing the system though they may appear 
to be in radically different form. The characteristic equations obtain- 
able by these methods for a given system should be identical. The 
technique discussed under the two separate headings either have been 
developed through the approach within which they are included or they 
are most frequently applied in that type of approach to the problem. 
Actually none of the techniques discussed is restricted to application 
by one approach or the other. 

For systems having complicated feedback hookups between the compo- 
nents Identified as the control unit and the aircraft, the resulting 
interaction of these components makes it impossible to treat them as 
cascaded elements. In such cases the combination of the autopilot equa- 
tion with the equations of motion is complicated and results in a higher 
order system of equations than that obtained by the simple substitution 
permissible with cascaded elements. The transient-response analyses -for 
such systems would probably be so Involved a,s to make the procediu’es 
impracticable except where an analogue computer were available. The 
frequency-response techniques, on the other hand, can handle the feedback 
complications rather neatly. (See references 3 and 22.) 

Autopilot as an added term in the equations of motion .— The types 
of analyses relating to this approach are essentially closed— loop analy- 
ses; that is, the motion of the aircraft-autopilot combination is 
analyzed with the system completely assembled and operative. Open— loop 
analyses, where the chain of operations is broken at some point in the 
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loop such as the output feedhack to the sensing device, are used in the 
servomechanism type of analysis. The open-loop analyses lend themselves 
more readily to the synthesis of automatic control systems than do the 
closed— loop type. (See reference 34, page 518.) 

The Routh or Hurwltz criterion ;— The simplest and usually most 
rapid test of the stability of a dynamic system is to apply the 
Eouth criterion. (See reference 37-) This criterion was put in 
the form of some easily written determinants by Hxarwitz. ' It 
determines whether or not all of the roots of the characteristic 
eq.uation have negative real parts, thus indicating whether or not 
subsidence of aperiodic modes and damping of the oscillatory 
modes exist. 

The first step in the application of this technique is to 
write down the characteristic equation in the following form: 


a^ + a^ D°-“^ + 


+ an = 0 


( 66 ) 


. t , an are real constants and a^ > 0. A neces— 


where a^, a^^, 

sary condition for all the roots of the above polynomial to have 
negative real parts is that all of the coefficients be positive. 


According to the Hurwltz criterion, the following determinants 
must be positive for a stable system: 
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J 

For a 'System of nth order there are n of these determinants . 

A derivation of- the Hurwitz criterion is given in reference 38 . 
In this reference, the criterion is stated in the following form: 
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"The roots of a polynomial with real coefficients 
have negative real parts if and only if the above 
determinants are positive, provided the coeffi- 
cient ao is positive." 

Eouth, in reference 37^ derived this same criterion. The 
derivation used hy Eouth is based on the theorem of Sturm. (See 
reference 25.) It is somewhat easier to use the Hurwitz form, 
since it can be written conveniently in terms of the deter- 
minants. These criteria are applicable to polynomials of any 
order . 


In most cases it is not necessary to compute all n of the 
Hurwitz determinants. Certain of them may be red\xndant. Con- 
.sider, for instance, the case of a fourth-order polynomial. 

The Hurwitz determinants are 


(ai ag - ao ag) > 0 
(a^ag - a^ ag) - a^>0 


(ay 

(b)' 


( 68 ) 


Since all of the coefficients must be positive, it is clear 
that condition (68a) can be derived from condition (68b). 

The determinant (68a) is therefore superfluous. In the case 
of high-order polynomials, much labor can be saved if super- 
fluous determinants can be found and eliminated. 


Numerical evaluation of roots ;— In general, the charac- 
teristic eqviations of linear systems are polynomial eqiiations 
that can be solved by well-known methods. No details of such 
methods will be given here , but the Graeffe root squaring 
method (reference 24) and an iteration process which can be 
carried out to obtain any desired accuracy (reference 3> PP- 
89-91) are recommended. 

In the case, of the transcendental characteristic eq\aa— 
tion, such as arises through the \ise of the. lag operator, a 
graphical method of solution is recommended. , Such a solution 
is given in reference 23. 

Stability boundary charts :— Basically this method is an 
extension of the Eouth or Hurwitz criterion method. If the 
influence of two variables on the stability of the motion is 
of interest, a plot can be, made showing stability boundaries 
with these two variables as coordinates. An example describ- 
ing the construction of such a chart is given in appendix A. 
The criterion for establishing the stability boundaries is the 
vanishing of the real parts of the roots of the characteristic 
equation. 
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Many different forms of such charts are employed in various 
fields of mechanics for the purpose of dynamic stability analysis. 
(See references 6, 26, 39# ^0# ^1# and 42.) For instance, if the 
values of the pefiod and damping of the motion are of more interest 
than the absolute stability boundaries, lines of constant period 
and damping may be substituted for the stability boundaries or 
even superimposed on them if desired. The superposition of such 
lines and boundaries could be extended to cover even more criteria, 
but the chart would soon become a confusing network of lines. 

Presumably any two independent variables of interest, usually 
relating to the physical properties of the system, such as stability 
derivatives, control gearings, moments of inertia, or a coefficient 
representing an algebraic combination of such parameters, can be 
used as the two coordinates. The stability characteristics for • 
which the boundaries or lines of constant magnitude are constructed 
are usually related to the output or behavior of the system. These 
characteristics axe often in the form of stability boundaries, 
damping and period constants, regions of optimum performance, etc. 

Examples of such charts and methods of construction are given 
in references 6, 26, 38# 39# 40, 43, 44, and 45. 

The theory underlying the construction of the stability 
boundary chart in particular, such as those constructed in 
appendix A, will now be discussed. 

The characteristic equation of the system glveS' the roots 
for all the possible modes of motion. Eouth (reference 37) and 
Hurwltz (reference 25) have determined that if the coefficient 
of the constant term of the characteristic equation is positive 
and if certain inequalities, developed by Eouth as discriminants 
and also obtainable from the Hurwitz determinants, axe satisfied, 
the real parts of all the roots will be negative and all modes 
of motion will be stable. (For proof of this see reference 25.) 

By equating this coefficient of the constant term to zero and 
by considering the inequalities to be equalities it is possible 
to establish the boundaries of a stable region in terms of any 
two independent variables involved in the coefficients of the 
characteristic equation. For any characteristic equation of 
polynomial form (for the case of transcendental characteristic 
equations, see reference 23) the number of determinants that 
exist is equal to n, the order of the equation. In general, 
then, it would seem that n+1 functions must be plotted to deter- 
mine the region of stability. Actually only two functions must 
be plotted to determine the stability boundary. The locus of 
an=0, where , a^ is the coefficient of the constant term, and 
the locus of the a— 1 determinant give the complete boundary. 



3k 


mCA TN 2275 


The locus of a.^=0 forms the houiKiary between stable and 
unstable aperiodic modes. This statement is based on the fact 
that in proceeding from a stable to an mstable condition the 
coefficient a^j which is eq.ua! to the product of all the roots, 
is the first to change the sign if one real root does but it is 
not affected by changes in sign of the real parts of the complex 
roots . 

The locus of the n— 1 determinant fonns the boundary between 
the stable and \mstable oscillatory modes. This can be shown^ 
by talcing lo), the value of the root to the characteristic equa- 
tion which would yield neutral oscillatory stability, and sub- 
stituting it for D in the characteristic equation. Then separat- 
ing the real and Imaginary parts and applying Sylvester's 
elimination theorem (reference 25) co can be eliminated from the 
separated equations leaving a determinant which is precisely the 
n— 1 Hvirwltz determinant . 

A variation of this method of determining the oscillatory 
boundary has been developed which uses numerical values of o) 
in the characteristic equation rather than the locus of the 
n-1 Hurwitz determinant. This numerical process seems to be 
somewhat simpler and provides, additionally, the values of the 
frequencies at various points along the boundary. Details of 
the construction of such a boundary are given in references 38 
and 46. The procedure, generally speaking, is to substitute 
the root ioj into the characteristic equation, set the real and 
imaginary parts equal to zero, and solve the resulting equations 
simultaneously (in terms of the two chosen Independent variables) 
for a sufficient range of values of oj to define the oscillatory 
boundary. Application of this method for obtaining lines of con^ 
stant damping and frequency is very easily carried out. Such 
procedures using the numerical value of cd are described in 
references 23 and 44. 

Hence, to summarize the technique of plotting the stability 
boundaries, it can be said that the oscillatory mode part of the 
stability boundaries can be determined analytically from the 
nr-1 Hiirwitz determinant or numerically using io) as a root of 
the characteristic equation. The portion of the boundary 
separating the regions of stable and unstable aperiodic modes 
is the locus of the coefficient of the zero-order term. 

Transient response analysis ;— If more information is desired 
than can- be obtained from application of the first three tech- 
niques discussed, the transient response, that is, the time 
history of the motion following a disturbance, can be investi- 
gated. Various kinds of step and pulse inputs are ordinarily 
2 

Proof communicated to the authors by Mr. Harry Greenberg. 
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applied in the calculations or experimentation used to obtain the 
transient response. The unit step function is perhaps the most 
commonly employed input. It is discussed at greater length in 
the section on evaluation of techniques. The transient response 
can be obtained analytically in at least four ways; (l) by the 
classical methods for solution of simultaneous differential eqiia— 
tions, (2) by the operational method of solution, (3) step-by-step 
calculations, and (k) by the use of an analogue computer. For the 
first two of these methods the roots must be known. Expres- 
sions for the roots can easily be written down for most second- 
order systems. For higher-order systems it is usually more prac- 
tical to solve for the roots nvimerically, which necessitates 
application of one of the techniques previously discussed. The 
roots need not be obtained in the application of the other two 
methods. They can be estimated roughly, however, from the result- 
ing curves showing the time histories of the motion. If the 
frequency— response data are available the transient response can, 
of course, be obtained by Fourier synthesis. 

The usual procedure followed in applying the first two 
methods is to formulate the characteristic equation, determine 
its roots, formulate the complementary and partlcul^ solutions, 
and insert the known boundary or initial conditions in order 
to evaluate the undetermined coefficients. Generally speaking, 
the character. of the transient response is determined by the 
characteristics of the servomechanism itself and is independent 
of the form of input function. Thus the particular solutions 
are used solely to determine the steady— state error, even though 
varioxxs transients result from various input functions, and the 
transient response is obtained solely from the complementary 
solution. 

The application of analogue computers to obtain the time-liistory 
solutions is a very powerful technique. A discussion of the details 
of using computers is not within the scope of the report. However, 
it is obvious that in order to carry out the calculations on such 
machines the numerical values of the constant coefficients and of 
the constants related to the initial and boundary conditions must 
be known. These analogue computers can be applied to solve linear 
systems of any order within the limits imposed by the size of the 
machine. Investigation of nonlinear systems is also possible. 

Examples of transient solutions using the airplane stability 
equations, with the autopilot included or at least approximated 
by a lag factor, can be found in references 23, 36, and 37. 

For fourth-order systems representing an uncontrolled aircraft, 
references ^7 and give transient response solutions. In 
general, however, the automatically controlled systems of 
fourth or higher order are more readily analyzed using the 
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frequency response techniques to be described later. Consequently, 
most of the solutions for such systems appear* in the frequency- 
response form. 

Airplane response as a transfer function in the control, loop .— The 
techniques to be discussed in this approach will be divided into two 
groups: the first representing transient response analyses; the other 

representing frequency response analyses. As was the case for the pre- 
viously discussed techniques, those in the transient response group will 
be basically closed— loop analyses; whereas those in the frequency response 
group will be essentially open-loop analyses. 

The results of the open— loop techniques are always capable of being 
interpreted in terms of the closed— loop performance, otherwise there 
would be no point in. carrying out such analyses. The transient response 
is the most pertinent for Judging the results of the analysis and, con- 
sequently, the open-loop and closed— loop frequency response techniques 
must provide results that can be related to the transient response. 

Transient response :— As mentioned previously the initial 
equation or set of equations developed for the dynamic analyses 
of servomechanisms may have been obtained by either of two 
formally different approaches considered in this report, but 
the eqixations so derived will be independent of the approaches 
except possibly for their mathematical form. Thus it is reason- 
able that the discussion of the technique of transient response 
analysis under the aircraft equations of motion approach will, 
in general, be applicable to the servomechanisnt-loop approach 
also. 

Examples of the use of transient response techniques in con^ 
junction with this servomechanism-loop approach are of interest, 
however, and. are worthy of mention. This technique is applied 
extensively in reference I 3 to the second-order servomechanisms 
with various types of control and damping. Classical methods are 
used therein to solve the differential equations encountered 
because, as the authors stated, in contrast to the more advanced 
process of operational calculus these methods more readily 
afforded an understanding of the physical phenon^na involved. 

Examples of both the classical and operational solutions for 
second-order servomechanisms are given in reference 3 * Refer- 
ence 20 applies the operational methods and gives a ^standard 
form" criterion for use in analyzing and modifying the stability 
and performance of the system. These criteria are established 
for systems of the order of two to six and are based on an 
initial overshoot of from 10 to 20 percent. A table of charac- 
teristic equations with numerical values of the coefficients is 
given. The equations listed correspond to various percentages 
of overshoot and are intended for use as a basis of comparison. 
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that is, a standard. The method, however. Is not necessarily 
restricted to those standard forms chosen in reference 20. If 
different performance characteristics are required the standard 
forms can he, suitably altered. A third-order system is analyzed 
by classical methods in reference 4l and by operational transient- 
response methods in reference 26. 

With the possible exception of the "standard form" technique, 
however, the fourth and higher order systems probably. can be more 
readily analyzed by the frequency— response techniques to be dis— 
ctissed. 

Frequency response :— Recently the method of frequency— response 
analysis has become a very popular method for studying servo- 
mechanisms. Basically this method is a study of the effects of 
a sinusoidal input of various frequencies on the amplitude ratio 
and phase characteristics between the input and the response, 
or output, which is sinusoidal also. The method is applicable 
to both experimental and analytical techniques. This versa- 
tility, plus the fact that it is a relatively fast method for 
purposes of synthesis, if not analysis, accounts in part for 
its popularity, 

A good introduction to this method is given in reference 3, 
pages 92—96, which includes discussions of certain general 
relationships between transient and frequency behavior and a 
treatment of the preparation of system data into a form con- 
venient for analysis and synthesis by frequency-response 
techniques. A short discussion of the various techniques is 
given in the following paragraphs . 

.(1) Resonance plots 

N 

(a) General applications 

If the sinusoidal input is written as 

®i ^ ' (69) 

the resulting output .will be of the form 

0 = Rg cos (cot + e) (70) 

where e is the phase angle between the output and input. These 
expressions may be considered as the real parts of 

01 = R^e^^^ 


and 


(71) 
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(72) 


by virtue of the well-known expansion. of a complex exponential 
function into sine and cosine components 

= cos ojt + i sin ost (73) 


With the input and output expressions written in this exponen-^ 
tlal notation it is easy to obtain a form for the closed loop 
transfer function of the system. Thus 


_0_ ^ Eg ^ ) 

®i " E^ 


(7^) 


or 


® 1^2 ie 


(75) 


Eg 

where — is the amplitude ratio. 


Plots of the amplitude ratio of the output to input 



9 _ 

9i 


and of the phase angle e against frequency u 


give the basic closed— loop form for the frequency— response 
characteristics of a system and are often called resonance 
plots. Some typical resonance cvirves for servomechanisms 
are given below: 




The frequency— response characteristics are most readily 
related to the transient response when presented in this form. 

It is not convenient to determine quantitative stability data 
such as critical control gearings, etc., however, from these 
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closed-loop resonance plots. It is possible, though, to deter- 
mine the existence of a steady— state error from the shape of 
the amplitude-ratio curve at 00 = 0 . For a displacement type 
of input such as a step function there is no, steady error if 

the frequency-response ratio of the amplitudes — (ioo) is 

1.0 and tThe phase angle is zero a*b (o = 0. Foi* zei*o steady— state 
error with constant speed or constant acceleration inputs^ there 
is the further rei^uireiiient that the first and second derivatives, 
respectively, of the amplitude-ratio curve and the phase--angle 
curve must he zero at zero frecjuency. 

Generally speaking, the evaluation of the behavior of 
closed-loop systems by frequency-^response methods starts with 
the determination of the open-loop response (discussed later) 
which is then related to the closed— loop resonance plots and 
to the transient response through the implicit correlations 
that exist between these techniques. 

The relationship between the sinusoidal behavior as 
indicated on the resonance plots and transient behavior 5s 
discussed at some length in reference 3. (See pages 104 - 106 .) 

To summarize -these relationships briefly, consider the follow- 
ing figure which has typical curves for the closed-loop 
transient response and the closed— loop frequency response for 
several values of the damping ratio t . These curves are 
representative of a second— order servomechanism,; but they are 
sufficiently general for illustrative purposes. 



The peak amplitudes on the frequency— response curves indicate 
resonance and the frequency at resonance for any given damping 
ratio ^ is a fairly good qxiantitative indication of the 
frequency of the transient oscillation. Also the magnitude of 
the peaks C€in be taken as a good indication of the relative 
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clamping of the transient response. For clamping ratios of approx- 
imately 0.6 or less the peak overshoots for corresponding curves 
can he fairly well correlated quantitatively. The speed of 
response to transient signals can be predicted from the frequency- 
response curves since a high resonant frequency would indicate a 
high natural frequency and, consequently, a high speed of response 
to an input signal. Thus the frequency-response ciorve for 

0.4 has its peak at the highest value of (n/oi^ of the curves 
shown and on the transient-response plot this value of ^ is 
seen to give the quickest response. The transient curve for 
C = 0.4 also shows the greatest overshoot as could he antici- 
pated from the relatively large magnitude of the peak on the 
C= 0.4 curve of the frequency response. For systems of order 
higher than the second it is possible that a high peak on the 
frequency-response curve will indicate a small oscillation in 
the transient-response curve before reaching 9 j|^ rather than 
an overshoot of It is usually possible to predict an 

overdamped transient response if the curve for the frequency 
response has no peak. For the cases illustrated, however, the 
curves for values of i between O. 7 O 7 and 1.0, show no peak on 
the frequency-response curves but the transient c\arves show a 
slightly underdamped, rather than overdamped, response. 

The phase shift between 9 . and 0j^ is also an important 
characteristic and the magnitudes of the phase shifts that are 
shown on the frequency— response curves can be used to explain 
Certain characteristics of the transient response. Limits are 
often set on these magnitudes in order to provide for suitable 
transient— response characteristics. 

( b ) Greenberg * s Method 

By plotting the resonance curves for the aircraft and 
autopilot separately, Greenberg (reference 7) able to find 
the critical control gearing for a system having all the other 
parameters fixed. The critical control gearing k^j. is 
specified as the value of the gearing between the autopilot and 
the aircraft control surface that provides neutral oscillatory 
stability (hunting oscillations) at the critical frequency. 

The critical frequency is the frequency at which the total phase 
lag of the aircraft-autopilot combination is either 0 ° or l 80 ° 
(depending on the definitions of the signs of the phase angles). 
The procedure for finding the critical control gearing corre- 
sponds to the method of finding the gain margin of the system 
which will be described later. 

The same plotting technique is used in reference 21 by 
Jones and Sternfleld, but the stability condition applied in 
this case was, essentially, that the total phase lag shall be 
less than l 80 ° when the product of the amplitudes Is unity. 
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This procedure corresponds to the method of determining the 
phase margin for the system as described later. 


(2) Polar plots .— The use of polar plots showing the magni- 
tudes and phase angles of the frequency— response characteristics 
by curves tracing the vector loci are possibly the most fre- 
quently used type of diagrams in automatic-control— system 
analysis at the present time. These loci show the variations 

9 


of the vector qmntities such as 


9 ^ 


0 E 

= or — over the frequency 

ht & 


range 0 < aj<<» as continuous curves plotted in the complex 
plane. Of the various vector quantities that could be plotted 
the vector locixs of more commonly known as the transfer 

® A 

function locus, is the most useful. The quantity ^ represents 

E 

the open-loop transfer fimction of the system. 


(a) Open— loop transfer locus plots (Nyquist diagrams) 

Q 

The open— loop transfer locus ^ is related to the over- 
all frequency— response characteristic of the closed— loop' system 

— in the following manner: 

®i 

9 


9/E 

1 + 9 /E 


(76) 


The open-loop response is often preferred to the closed- 
loop response in both analysis and synthesis because it is more 
easily obtained and is more useful in studying the effects of 
system parameter changes on the closed-loop stability character- 
istics. The closed-loop response can be obtained readily from 
the opei>-loop response ciirves by a simple graphical relationship 
that is demonstrated on the sketch of the open-loop transfer 
locus plot. 

In using this type of plot 
to determine the behavior of a 
system it is convenient to consider 
three separate regions of the locus 
and inspect them individually. 

These regions are the low-€requehcy 
region, the region near the point 
— 1 + Oi, and the high-frequency 
region near the origin. 

From the low-frequency 
region it can be determined whether 
the system is a zero— position error, 
a zero— velocity error, or a 
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zero-acceleration-error servomechanism according to the axis 
the locus approaches at zero frequency. These classifications 
are related to the mamher of integrations in the loop and 
their significance is discussed at some length in references 3 
(pages 166— 168) and 20. 

The region near the point —1 + Oi relates to the sta- 
bility of the system and it is by far the most important region. 
The high-frequency region may possibly yield information 
regarding the over-all complexity of the physical system. 

From the —1 + Oi region it can be determined whether 
or not the system is stable by use of . the Nyquist criterion^ what 
the gain and phase margins are^ and some indication of the magni- 
tude of the peak response^ which relates to the maximum overshoot^ 
and of the resonant frequency. 

The Nyquist criterion can be best explained by relating 
it to the roots of the characteristic equation 1 + 9/E = 0. 

It was stated previously that if any of the real parts of the 
roots of this equation were positive the system was unstable. 

This is equivalent to saying that all the roots^ when plotted 
on the complex p plane, must lie to the left of the imaginary 
axis for stability. Foots lying on the axis give neutral sta- 
bility and roots lying to the right of the axis indicate diver- 
gence of the motion. The Wyquist criterion, however, relates to 
the frequency— response plot of the transfer locus on the 0/E 
plane* The* relationship between the transfer locus' plot and the 
location of the roots on the p plane is Involved with the mapping 
of regions from one plane to the other. It can be shown that 
the imaginary axis of the p plane, defining the boundary between 
the stable- and unstable roots, maps as the transfer locus on the 
0/E plane. Furthermore, as .shown in references 3, 11 , I3, and 
49, if the transfer locus encloses the point —1 + Oi there are 
roots on the right-hand side of the imaginary axis of the p 
plane and the system is unstable. If the transfer locus passes 
through the point —1 + Oi neutral stability is indicated. 

In the manner just stated, this criterion is rather vague and 
applies only to simple systems having no multiple loops. A more 
specific "rule of thum'tJ' for the stability of single— loop systems 
is given in reference 13 - The rule says, in effect, that, if a 
line originally coinciding with the imaginary axis on the 0/E 
plane can be fitted to the transfer locus through a clockwise 
rotation and a suitable distortion without sweeping over the 
— 1 + Oi point, the’ system is stable. Another quick test for 
single-loop systems is to trace along the locus in the direction 
of increasing frequency and if the —1 + Oi ^point lies to the 
left of the curve as the tracing proceeds to the origin the system 
is stable. 
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These simple tests may need modification when the system 
•is complicated hy multiple loops. In general, the Nyquist 
criterion requires that for stable systems, that is, systems 
having no positive real parts in the roots of the characteristic 
equation, the number of times that the locus encircles the -1 + Oi 
point shall be equal to the number of poles of the transfer 
function. Such poles may exist in a complicated system inasmuch 
as one of the subsidiary loops may be unstable by itself, while 
the closed-loop system is stable. A proof of the Nyquist 
criterion covering multiple-loop systems is given in reference 11 
and references 3 and k 9 give an extensive treatment of the subject. 

The phase and gain margins give some idea of the degree 
of stability, although the degree of stability can be found only 
by knowing the values of the roots of the characteristic equation. 
The gain margin, designated as , Gin in the following figure, is 
defined as the reciprocal of the magnitude of the transfer- 
function vector at the phase angle of -180°. Defined in this 
manner the gain margin can be interpreted easily as decibels 
below the zero decibel level on the logarithmic plots to be dis- 
cussed. Another definition of gain margin, namely, the differ- 
ence between the magnitude of the transfer function and the 
-1 + Oi point at a phase angle of - 180 °, has also been used, but 
it is not a practical concept. 

In cases where the transfer locus crosses the negative 
real axis more than once there can be as many values of gain 
margin as there are crossings. This situation occurs in the 
illustrative example presented in appendix A. The phase margin, 
designated as in the following figure, represents the 

amount by which the phase of the 
transfer function falls short of 
being — l80° when the amplitude of 
the transfer function vector is 1. 

, An indication of the 
magnitude of the peak response and 
of the resonant frequency can be 
obtained from the transfer locus 
plots by superimposing contoiur lines 
of constant magnitude and constant 
phase for the Q/q^ -vector. 

Examples of such contour lines 
designated by M and W, respec- 
tively, are shown on the following 
‘chart: 
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It Is otvious that by use of these M-N contours it is possible 
to mark off zones on the plot through which the transfer locus 
must pass in order that the magnitude and phase lag of d/d\ 
shall fall within prescribed bounds. Furthermore, the M 
circle which is tangent to the transfer locus gives the peak 
magnitude of the response and the frequency at which 

this occurs will be the resonant frequency. 


The most obvious application of these contours is in 
the reshaping of the locus to obtain better response character- 
istics-. By vectorial addition on the locus plot the effect of 
inserting a new component in the loop can be determined fairly 
easily. The effects of increasing the gain of the system are 
qviite readily observed on the polar plots. 

(b) The inverse open— loop transfer locus and others. 

Other polar plots of various vector quantities relating 
to the open— loop or closed— loop response are sometimes employed. 
The inverse transfer locus E/0 and the closed— loop— system 
response locus 9 /Q^ are seen fairly often. 


The inverse transfer locus 
shown below has one advantage in that 
the 9^/9 response can be determined 
by simple visual inspection rather 
than by graphical construction. (See 
references 22, and 32.) The use 
of one or the other of these two 
locus plots, however, is largely a 
matter of experience and personal 
preference. 


CO increasing 


Imaginary 

axis 



( 3 ) Logarithmic plots .— The tech- 
nique of employing logarithmic plots 
of the amplitude ratio and phase 
shift against frequency is useful 


e/ 0 plane 
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primarily in synthesis, that is, in the prediction of the effects 
of adding corrective networks to the system. In this method the 
phase angle and the log of the amplitude ratio of the transfer 
fimction (usually given in decibels, db, which is defined as 20 
times the common logarithm of the amplitude ratio) are plotted 
against the log of the frequency as shown in the following sketch: 



Log d) Log 0) 

The use of logarithms makes it possible to add ordinates 
directly rather than vectorially when altering the response 
curve to take into account the addition of some component to 
the loop. Furthermore, the log moduli of first-order components 
of the loop plot either as straight lines or approximately as 
straight-line segments. Second-order components can be handled 
conveniently using previously prepared families of logarithmic 
curves with the damping factor as the parameter. Thus servo- 
mechanisms made up of first — or second— order components in 
cascade can be rather quickly plotted on the log modulus graphs. 
There is also a special way of handling the phase— angle plots 
that permits rapid manipulation of these curves. Very good 
descriptions of the logarithmic plot techniques are given in 
references 3 and 11. 


In order to relate the 
open— loop characteristics 
of the transfer function 
to the closed— loop charac- 
teristics a cross plot can 
be made with the magnitude 
of the transfer function 
in decibels plotted against 
the phase angle for various 
frequencies. On this type 
of plot the M and N contour 
lines can be superimposed 
to determine the closed— loop 
characteristics. A sketch 
of this type of plot is given. 



-180 0 

Pliase angle, deg 
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The point of zero decibels magnitude and — l 80 ° phase shift now 
becomes the instability point to be avoided. 

The. manner in which the gain must be changed to stay within 

a certain limit of [ ^ | at resonant frequency and the 

°i 

max 

effects of adding components and changing the transfer-function 
locus shape can be determined on these logarithmic-coordinate 
plots in much the same manner that it was on the polar plot. 

The gain margin and phase margin also can be determined from 
these plots. 

This method of analysis can be applied to systems having 
higher than second-order components as is demonstrated in 
appendix A; but the facility of the method is definitely reduced 
luiless the system can be broken down into first- and second-order 
elements. In the case of the longitudinal motion of an aircraft 
the fourth— order equation representing the aircraft can be 
approximately factored as shown in appendix A. It is relatively 
easy to plot the quadratic factors and the example in appendix A 
shows that the results for the approximate factorization are not 
appreciably different from the curves for the quartic equation. 


Two Degrees of Control 


As Indicated previously the airplane, considered kinematically, is 
a slightly more complex object to control than those considered in the 
ordinary treatises on servomechanisms. Having six degrees of freedom, 
the airplane usually requires as many as three degrees of control to 
get satisfactory stability and performance. It is a well-established 
fact, however, that the motion of the conventional airplane can be 
divided into two categories, longitudinal and lateral, each with three 
degrees of freedom. The longitudinal motion, usmlly controlled by the 
elevator only, has been covered in the previously discussed analyses . 
pertaining to one degree of control. The lateral motion, on the other 
hand, is usvially controlled by some combination of the aileron and rudder 
deflections. In the lateral case it is theoretically possible to make 
the deflections of both the aileron and the rudder functions of the 
displacements and velocities of roll and yaw, of the sideslip velocity, 
of the accelerations in roll and yaw, etc. Practical limitations, how- 
ever, have so far restricted the analysis of lateral controls to two 
cases: (1) the so-called "cross— coupled" control (see reference 6) 

where the aileron and the rudder deflections are proportional to both 
displacement in bank and displacement in yaw (i.e., azimuth); and 
(2) the so-called "simple" control where the aileron deflection is pro- 
portional to displacement in bank and the rudder deflection is propor- 
tional to displacement in yaw (azimuth). Even in the last and most 
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simplified case, however, the analysis of the stahility of the over-all 
lateral motion is greatly complicated due to the Inherent coupling of 
the three individual lateral motions. The existing theoretical treai^- 
ments of the, two degrees of control arrangement are restricted to the 
latter case except for one "brief treatment of cross-coupled control in 
reference 6 . 


Autopilot as an added term in the eq,uations of m otion;— As in the 

case of one degree of control, this approach leads to the closed-loop 
type of analysis. 

Imlay*s method .— In this method the autopilot is introduced 
through a forcing function added to the equations of motion of the 
airplane. In the particular application of this method "by Imlay, 
reference 6 , an approximate delay time expression is used to simulate 
the action of the autopilot. This artifice is not essential to the 
method and with the now existing Imowledge of servo systems it is 
possible to write expressions that will give a good linear representa- 
tion of their action. 

Imlay considered the "cross-coupled" control "briefly. In spite 
of the complexity of the expressions resulting from this type of 
control he was able to determine the control gearings to give the 
special case of equal roots; This solution corresponds to motion 
with no oscillatory components and having all the aperiodic modes 
eqmlly damped. It was a .satisfactory solution at the lift coeffi- 
cient Investigated but the same control gearings caused instability 
at higher lift coefficients. Further study of this control was 
dropped in favor of a more complete investigation of the simple 
control. 


Simple control was investigated by the technique of developing 
stability boundary charts. The two variable parameters upon which 

the chart was based were the two control gearings, for the 


r\xdder, and 


a6, 


for the aileron. This technique involves the 


establishment of criteria in the form of algebraic inequalities by 
the Eouth or Hurwitz methods and the subsequent study of these 
inequalities in terms of the two variable parameters in order to be 
able to chart the stability boundaries. For a system represented by 
a fourth degree equation, such as the airplane uncontrolled, this 
procedure is fairly long and tedious. For systems involving higher 
order equations, such as an aircraft-autopilot combination yielding 
a sixth-order equation, the procedure may become so laborious as. to 
be prohibitive. 


Airplane response as a transfer function in the cont rol loop:— In 
addition to Greenberg’s method discussed below, the open— loop polar- 
plot technique (Nyqulst diagrams) could undoubtedly be applied to the 
analysis of aircraft-autopilot systems having two degrees of control. 
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A system using this type of control falls into the category of multiple- 
loop systems which the Nyquist diagram technique can handle, in prirt- 
ciple, no matter how complicated the system might he. The details of 
the application of this technique to the particular case of an aircraft- 
autopilot comhination having two degrees of control are not available, 
however, at the present time. 

Siii 5 )lified methods of handling two degrees of control are desirable 
and perhaps could be developed by neglecting the coupling of the control 
effects as a first approximation or by the use of some similar artifice. 
There is also room for the development and presentation of more powerfiil 
analytical techniques, especially if the time required for their appli- 
cation can be held within reasonable bomds. It is possible, however, 
that the use of a simulator, such as an analogue computer, will always 
provide the best combination of simplicity and versatility for the 
analysis of systems having two degrees of control. 

Greenberg's method.— This method (reference 50) is application 

of the frequency— response type of analysis to the calculation of the 
stability of an airplane with two independent autopilots or, in 
Imlay's terminology, with "simple" control. It is an extension of 
Greenberg’s earlier paper, reference 7^ on the application of 
frequency— response methods to obtain the critical control gearing for 
one degree of control. 

The method makes use of a type of resonance plot similar to that 
found in reference 7. The calculations involved are, naturally, much 
longer and more tedious than those for one degree of control. Basic- 
ally the method gives only discrete points of the stability boundary 
using the control gearings as the variables, but it could be applied 
a number of times in order to get a fairly accurate layout of the 
boundary. Whether this would be more or less time consuming than 
Imlay’s method is questionable. 

EVALUATION OF TECHNIQUES OF AMLYSIS 

Although no specific performance or stability standards are to be 
set up in this report, an attempt will be made in the present section to 
evaluate the techniques of analysis on the basis of their usefulness in 
providing information about the parameters affecting performance and 
stability. In the discussion of fimdamental analytical concepts the four 
parameters chosen for consideration were (1) speed of response; 

(2) steady state error; (3)' time to damp; and (4) amount of initial 
overshoot. It was pointed out that these items could be reasonably, 
though not absolutely, paired off under the headings of control and 
damping parameters. Furthermore, the damping and control can be ration- 
ally paired off with stability and performance, respectively, although, 
here again, the division is not absolute. Thus the following diagram- 
matical outline is developed relating the four parameters to the behavior 
of the system in terms of performance and stability: 
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I. Routh’s Discriminant 

This technique merely provides the answer as to whether stability 
does or does not exist. Thus it is only useful for obtaining qualita- 
tively information relating to parameter 3* 

II. Numerical Evaluation of Roots • 


The numerical evaluation of the roots provides information regard- 
ing the time to damp (parameter 3). Furthermore from the complex roots, 
the natural frequency of the oscillatory modes can be calculated and 
used to estimate roughly the speed of response (parameter 1) . 

III. Stability Boundary Charts, 

Basically the use of stability boundary charts permits combinations 
of various values of two Independent variables to be charted info regions 
of stability or instability. In this respect it is slightly more flex- 
ible than technique I but gives essentially the s^e information. It is 
possible, however, to superimpose a lattice of lines of constant damping 
and period on such a chart. With this supplementary feature, it is 
possible to use the charts to determine quantitatively the time to damp 
and to estimate the speed of response as can be done using technique II. 

IV. and V. Transient Response 

From a graph showing the transient response, information concerning 
all four parameters can be obtained. The determination of the transient 
response by theoretical procedures, however, such as described in the 
discussion of the techniques of analysis, is usually restricted to cer- 
tain standard types of input signal. The information obtained, conse- 
quently, might seem to be relevant only for operating conditions u nd er 
which such signals or similar disturbances occurred. Experience has 
shown, however, that this procedure is a satisfactory one and can be used 
rationally for analysis and synthesis of servomechanisms. , The response 
to a step input, for instance, can be used to evaluate the response to 
random signals or disturbances, or it can be applied, using mathematical 
techniques, to determine the system response to any arbitrary command. 

If the transient response of the servomechanism can be measured 
under actual operating conditions or by simulators, a rigorous analysis 
of the system behavior can be made. 

In addition to analyzing the transient response to a prescribed 
input or to a random input corresponding to actual operating conditions, 
there is the possibility of using the frequency-response analyses to be 
evaluated subsequently. In regard to the relative merits of the two 
types of analyses, reference 3 presents the following views, ".... that 
dynamic stability specifications are more' easily based on transients 
following a^ sudden disturbance than on the behavior following sinusoidal. 
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disturbance of varying frequency. However, the design of conventional 
systems to fulfill the majority of performance specifications is more 
easily accomplished by studying the response to a sinusoidal disturbance 
of varying frequency.” , 

A particular technique involving the transient response and yet not 
requiring the time history of the behavior of the mechanism to be 
ascertained is worthy of mention. Whiteley in reference 20 uses "stand- 
ard forms" of -the characteristic ^equations to specify the suitability of 
the system behavior. The standard forms are based on a given percentage 
of initial overshoot, a nonoscillatory response, and a satisfactory 
speed of response. Thus by comparing the numerical constants of the 
characteristic equation for any servomechanism with one of Whiteley’s 
standard forma, some idea of changes in ‘the constants, and thus in the 
system parameters required to match the performance of the servo- 
mechanism with that corresponding to the standard forms, can be obtained. 

VI . Eesonance Plots 

The resonance plots are usually employed to show the closed— loop 
frequency— response characteristics of a system. Information about the 
parameters under consideration, however, can be obtained only Indirectly 
from these plots. The speed of response (parameter 1), for instance, 
can be estimated from the location of the resonant peak through the 
relationship between resonant frequency and the natural frequency which 
has a correspondence to, the speed of response. The existence of a 
steady-state error (parameter 2) is indicated in this type of plot. The 
time to damp (parameter 3) and the amount of. initial overshoot (para- 
meter k) can be related to the magnitude of the peak response as des- 
cribed in the discussion of this technique of analysis. 

Methods using open-loop resonance plots, as done in references 7 
and 21, can be applied to determine gain or phase margins or related , 
parameters such as the critical control gearing. 

VII. Polar Plots 


The Nyqulst diagram (open— loop polar plot) is by far the most 
commonly employed type of polar plot. It is fairly convenient to con- 
struct and easily read by experienced analysts. Here again the informa- 
tion regarding the four behavior parameters is obtained indirectly 
through the correspondence of certain characteristics of the transfer 
locus with the transient behavior. The time to damp in terms of a value 
for a damping factor is not obtainable. However, quantitative informa- 
tion regarding the stability can be. obtained in terms of gain or phase 
margins or critical control gearings. 

Contour lines of constant magnitude can be superimposed on the 
Hyquist diagrams to determine the control gearing required for any 
desired peak response of the system. The peak response can be related 
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to the damping and initial overshoot as previously mentioned. If the 
frequency spectrum is indicated on the locus, the resonant peak location 
will give some indication of the resonant frequency which relates to the 
natural frequency and speed of response. 

The type of steady-state error can be estimated from the position 
of the asymptote of the transfer locus (or the mathematical form of the 
transfer function) at low frequencies. (See reference 1.) The type of 
steady— state error is often used as a means of classifying the servo- 
mechanism as a zero— position error, zero— velocity error, or a zero- 
acceleration— error system. (See references 3 and 20.) If a steady— state 
error exists, its magnitude can be determined best from the differential 
eqiiation for the system. 

7111 . Logarithmic Plots 


Essentially, everything that was said with regard to the information 
obtainable from the polar— plot technique applies to this technique as 
well. By enabling operations that required multiplicjation or division 
in the polar— plot technique to be performed by addition or subtraction 
and by breaking the transfer loci down into combinations of first— and 
second— order components that can be easily remembered and readily 
applied, the logarithmic— plot technique was developed as the most rapid 
procedure for analysis and synthesis. For synthesis there cannot be much 
argument about the improvement in facility gained by the use of logarith- 
mic coordinates. From the analysis standpoint, the advantage of using 
logarithmic charts is questionable. The open-loop transfer loci on the 
log amplitude versus log frequency and phase angle versus log frequency 
charts give data which, unlike the polar plots, cannot be readily inter- 
preted in terms of the closed— loop response for purposes of analysis. 
Therefore, another chart having an ordinate scale of the log modulus and 
an abcissa scale of the phase angle must be prepared. The transfer locus 
is then plotted and contours of constant amplitude and phase angle of 
the closed— loop response can be superimposed. From this type of log plot 
most of the stability and performance characteristics can be determined. 

IX. lTTilwy*a Method 

This technique (reference 6) for analyzing two-control systems is 
a combination of the numerical root solving technique and the stability 
boundary chart technique previously discussed. The independent variables 
selected for the latter are , usmlly the control gearings. 

X. Greenberg's Method 


This technique (reference 50 ) is an extension of the open— loop 
resonance plot technique to the two-control case. It is used mainly for 
determining critical control gearings and it is the only frequency- 
response technique reported that, at present, has been applied to an 
aircraft having two degrees of control. 
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CONCLUDING REMARKS 


The commonly employed techniques of stability analysis for automat- 
ically controlled aircraft have been discussed and classified in this 
report. The survey was limited to continuous control and to linear 
closed— loop systems. A discussion of the components and concepts asso- 
ciated with controlled aircraft and similar servomechanisms was included 
in order that the descriptions of techniques could be made fairly con— 
cise. An attempt was made to evaluate the techniques on the basis of the 
kind and amount of information desirable. Furthermore, an illustrative 
exEimple is included in an appendix that presents calculations for a 

typical case using all of the techniques discussed. 

« 

No attempt was made to establish satisfactory standards for the 
performance and stability of automatically controlled aircraft. Foirr 
performance and stability parameters were selected, however, to be used 
as a basis for comparison and evaluation of techniques. These parameters 
were: (1) speed of response, (2) steady— state error, (3) time to damp, 
and (4) amount of initial overshoot. 

The evaluation of the techniques can be sxanmarized briefly as 
follows: For a simple determination of system stability or instability,* 

application of the Routh discriminants (or Hurwitz determinants) is a 
sirfficient and relatively fast procedure. If knowledge of the rate of 
damping is desired in addition to the determination of stability, the 
roots of the characteristic equation can be calculated. When the design 
of the system is set except for the value of the control gearing to be 
used, application of the critical controL- gearing method or the polar 
plot of the open— loop transfer locus is recommended. If the system is 
likely to require more adjustment than the setting of the control- 
gearing value, then the polar plot is the more suitable of the two methods 
Use of the technique of plotting the transfer locus using logarithmic 
coordinates seems advisable only when a full-scale synthesis Job is 
required. Analyses sufficiently complicated to warrant application of 
an elaborate technique, however, are more simply and intuitively, handled 
by the polar plots. For the simultaneous investigation of two inde- 
pendent design variables with regard to system stability, the stability 
boundary charts are recommended. For the most elaborate and detailed 
analysis, that is, one that will provide information about each of the 
stability and- performance pareimeters considered, a transient response 
analysis of some sort is the best-. Transient-response techniques, 
however, are not as adaptable to synthesis as the frequency-response 
techniques . 

Many of these techniques are only applicable to cases where one 
degree of control is employed. Extension of these techniques to cover 
the situation where two degrees of control must be considered, such as 
in the case of the lateral stability of an aircraft, has been done in 
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only two generally available references, 6 and 50* A discussion of these 
techniques is included in this report. The possibility of using simula- 
tors for solving the two-degrees-of-control case is mentioned also. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 5, 1950. 
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APPENDIX A 
ILLUSTRATIVE EXAMPLE 


For the purpose of demonstrating several of the techniques of sta— 
hillty and control analysis described in the body of this report, calcu- 
lations have been made for the case of a typical high-speed fighter 
aircraft having an automatic pilot with proportional control based on 
the error in pitch angle E where E = 9^-0. The following diagram 
defines the variables involved and their positive directions; 



where 

•V . aircraft velocity 
u longitudinal, perturbation velocity 
ttg initial angle of attack 
a perturbation angle of attack 
9 angle of pitch ' 

9j^ input angle of pitch 
m mass of the airplane 
g acceleration due to gravity 

For this arrangement the longitudinal equations of motion of the airplane 
in level flight with controls fixed can be written (for a derivation of 
a similar set of equations , see reference 39) 
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or-2 T Dol+2 t D0+ ^ sin 


a„-S^\' = 




(■^To ) 

CL^^a+ (ci^-tCj^ sin aQ^0+ (2Cp^-2CT^ cos ao + 3^)^’ 

fl. ^ + 2^ 

^ T hxi^* ) 


0 (Al) 


T Du = 0 


^C] 

■§s:“ ■54 




(A2) 

(A3) 


where 

Ct 

^L 

Cd 

Cm 

K' 

u* 


thrust coefficient 
lift coefficient 
drag coefficient 
pitching moment coefficient 
nondimensional moment of inertia parameter 

nondimensional longitudinal perturbation velocity 


T aerodynamic time 

/ 

a, 6, etc. ^ ^ etc. 

dt' dt^ 

0 ^ 
dt^ 


p^o) 


me / 

(^) 


, seconds 


and subscript o denotes initial condition. 

For an altitude of 10,000 feet and a Mach number of O .85 the numerical 
values of the mass and aerodynamic constants are assumed to be 


( 003 oo - - - 5.325 (^ ) • -0.0165 

( _aci,^^CT^_ ^) = 5.39 ._o.o34 
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0.1553 



+ Cm 
■^0 


3ln a 




— 2Cm COS 
■^o 


= 0.0687 



0.553 


f i ^ 



2 


^ = -0.639 

000 


=^-0123 

^rn 


T = 0.967 




-0.6299 


Substituting these numerical values into equations (lA), (2A), and (3A) 
yields 

(-5.325) a - 2(0.967) Da + 2(0.967) D9 +(5.39)u* = 0 (A4) 

(0.155) a +(0.069)9 + (0.553) u' + 2 (0.967) D u* = 0 (A5) 

(-0.639) a + (-O.Ol65)Da + (-0.03^) D9 + (-0. 0123 ) (0.967)^0^9 + 

(0.1530) u' = 0 (A6) 


where D = — . 

dt 

It is also necessary to choose some form for the' autopilot equation. 
An electric motor having the usual inertia and viscous deimping charac- 
teristics and an internal feedback was selected. The expression for 
this autopilot is 

[ Jm D^ + (fm + S 3 ) D + St ] 6 = St k E (A7) 

where 

moment of inertia of motor 
damping coefficient of motor 
Sg torque— speed sensitivity of motor 

Sip torque— voltage sensitivity of motor 
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6 control deflection 

k control gearing 

With the numerical values assumed for these motor constants, the auto- 
pilot eq.uation becomes 

(0.058 D 2 + D + 37.5) 6 = 37-5 kE (A8) 

For this problem k was chosen as one of the independent variables to 
be investigated. The other variable chosen for the calculation of a 
stability boundary chart was 8Cjji/8a. This derivative, often referred 
to as the static stability factor, previously was set equal to -O.639. 
The loop diagram for the autopilot-aircraft combination is as follows: 



Having the equations repre- 
senting the aircraft and auto- 
pilot behavior, it is now neces- 
sary to determine the algebraic 
signs that will tend to make 
the combined action of the two 
mechanisms stable. According 
to standard aircraft nomencla- 
ture, the displacements shown 
on the sketch are positive. 

Thus a positive control deflec- 
tion 6 gives a negative angle 
of pitch 0 for the iisual 
aerodynamic forces and moments 

involved. Since the elevator deflection for the conventional aircraft 
gives merely a pitching moment about the center of gravity, only the 
third equation of motion will be affected by the inclixsion of the 



elevator action; With the elevator action 
function, equation (A3) becomes 


IT 


included as a forcing 
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^ a* ^ to. ^ ^ ^ |ffi u'= ^ ^ (A 9 ) 

or nmnerically 

(- 0 . 639 )a+(- 0 . 0165 )]>x+(- 0 . 034 )D 9 +(- 0 .0123 )( 0 . 967 )^I) 20 +( 0 . 1530 )u* = 

-6 (-0.6299) (AlO) 

As written, the autopilot eq.uation (a 8) gives a positive control 
deflection for a positive error since k is a positive quantity. In 
servomechanism nomenclature, the error is defined as 

E = - a (All) 

The following diagram represents the above relationship. Hence, to 

summarize, if a positive error 
gives a positive elevator deflec- 
tion via the autopilot and a 
positive elevator deflection gives 
a negative Increment to the pitch 
angle via the aerodynamics of the 
aircraft, the net result is an 
Increase in the positive error and 
the system so defined is static- 
ally unstable. Prom the aeronaii- 
tical engineer’s point of view, 
the easiest way to repair the 
system is to make certain that the 
autopilot gives a negative deflec- 
tion for a positive error, thus 
changing equation (A8) to read 

(0.058 D® + D + 37 - 5)5 = -37.5 k E (A12) 

Using equations (Al), (A 2 ), (A 9 ), and (A 12 ), the system is stable and 
completely represented in mathematical form. 

The results obtained by the applications of the various techniques 
demonstrated in the following sections of this appendix are summarized 
in table II. 


The Eouth or Hurwitz Criterion 


By replacing 8 in equation (A 9 ) with the equivalent expression 
obtained from equation (A 12 ), the equations of motion of the aircraft- 
autopilot become 
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cos clq — rDa+2 t D0+ ^-2Cj|^ -2C^ sin aQ — ^^^u'=0 (Ai) 


cos + 




u* + 


(A2) 


(f 5^ * ° 

E ^ 

.058 D^ + D + 37.5 / (AI3). 






Upon replacing the coefficients in the above equation by the munefical 
values, as was done previously, and setting k equal to unity, the 
equations become 


-5.325 a - 1,93^ DO-+ 1.93^ D0+ 5.39 u* = 0 


0.1553 a + 0.0687 ® + 0-553 u' + 1.934 Du* = 0 


(A14) 


(A15) 


-O. 639 cjr-O.Ol 65 Dc 30 -O .034 D 0 -O.O 1153 D^ 0 + 0.153 U* = ^ 37 - 5 ) ( , 0 - 6299 ) E ^ 

r\r- O .TS. -NP-.r- ' ' 


0.058 D^+D + 37.5 


As a further simplification, let 9^=0, so that E = -0. The 
characteristic equation of the above system is 


(-5-325 - 1-934D) 

0.1553 

(-0. 639-0. OI65D) 


1.934D 

0.0687 

f -0.034d-0.01153D^ - 

^ J 3 . -62 ^ 

^5 y 


0.0581)2 +D+37.5 


5-39 

(0.553 + 1-934D) 

0.153 


= 0 


(AI7) 


The expansion of this determinant results in the sixth— order polynomial: 

0.007OOD® + 0 .i 64 d^ + 5-59 D* + 39-8d® + 5170^ + 80ID + 219 = 0 (ai 8) 
It is seen that all of the coefficients of equation (AI8) are positive. 
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The independent Hiarwitz determinants for the sixth-order, polynomial 

eire 

0.I6I+ 39.800 

= 0.638 >0 (A19) 

0.007 5.590 

0.164 39.800 801.000 

0.007 5-590 517.000 = 12.41 > 0 (A20) 

0 0.164 39.800 

0.164 39.800 801.000 0 0 

0.007 5.590 517.000 219.000 0 

0 0.164 39.800 801.000 0 

0 0.007 5.590 517.000 219.000 

0 0 0.164 39.800 801.000 

Since the Hurwitz determinants are all greater than zero^ the 
system is stable for the given values of the parameters.' 

Numerical Evaluation of Roots 

The characteristic equation of the airplane-autopilot combination 
given by equation (AI8) is a sixth-order polynomial. The roots can be 
obtained by the method described in reference 3 and outlined in the 
following paragraphs . 

After division by the coefficient of D® the characteristic eqm- 
tlon becomes 

D® + 23.43 D® + 798.6 S'* + 5686 D® + 

73860 + 114400 D + 31290 = 0 (A22) 

The last three terms divided by the coefficient of D® yields 

+ 1.549 D + 0.4236 

as a first approximation to a quadratic factor. Division of the 


= 3,135,000>0 

(A21) 
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polynomial by this factor results In a remainder of 

3282 D + 3227 

and a second approximation of 

D2 + 1.699 D + 0.4723 

I 

obtained from the next to the last remainder. Division by the second 
approximate factor results in a remainder of 

105.0 D + 93.00 


Convergence is fairly rapid and only three more repetitions of the 
process are required to reduce the remainder to 

-39.00 D + 3.000 

This is small enough to be accepted as negligible. The corresponding 
factor is 

/ 

D^ + 1.700 D + 0.4737 


and the roots obtained from it are = -O.35I and ~ These 

real roots correspond to the damping constants of the two aperiodic modes 
shown in table I. , 

The remaining quart ic is 

D^ + 21.73 D® +. 761.2 D^ + 4382 D + 66050 


which can be factored by the same process to yield 

D^ + 4.097 D + 115.1 


and 


D 2 + 17.64 D + 574.0 


The roots obtained from these two quadratic factors are -2.05 ± i 10. 5 
and —8.82 ± i 22.3. Designating these as the roots for the long— and 
short— period oscillatory modes, the resulting damping constants and 
actual frequencies are 


( Coij,) = 2.05 

u 


“Z = 



10.5 radians /sec 
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( ^ = 8.82 


CDg = 03n >J 1- i = 22.3 radlans/sec 


which yields 


5 ^ = 0.195 

= 10,8 radiana/sec 

i g = 0.395 

cSfj = 24.3 radians/sec 

In this example the time to damp to half amplitude will be governed by 
the long— period root. Thus 

om , 0^ - 0 ^4 sec 

(EVj 2.05 


Cases for which the process is divergent can be handled by first 
transforming the equation by Graeffe’s root squaring method or by 
Horner’s root diminishing method. (See references 2 k and 25.) 


Stability Boundary Charts 


A stability boundary chart will be constructed using the control 

gearing k and the so-called "static— stability" derivative S'® the 

Independent variables. In determinant form, corresponding to equation 
(AI7), the characteristic equation for the aircraft-autopilot combina- 
tion is (with 9 ^ arbitrarily set equal to zero) 


(-5.325 - l. 93 ^D) 

0.1553 

( ^ - 0.0165d) 


1.93^ D 
0.0687 

-0.03l^D-0.0,1153D^ - 
23.62 k 
0.058 d^ +D+37 


5.39 

(0.553 + 1.93^D) 
0.153 




=0 (A23) 
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which reduces to the following polynomial form, 

aoD® + aiD® + agD** + agD® + a^D^ + a^D + ag = 0 


(A24) 


j'or t/h© C£l 36 cons id. ©1*6(1 "tli© co©f*f icion'ts of "bh© polynomial tscom©^ in 
terms of k and Cuj^ ^a symbolic notation for SCj^/Sa 


^0 = 

0.0000700 


ai = 

0.00164 



0.05230 - 

0.00576 Cin^ 

^3 " 

0.3340 - ( 

).1011 Cjn 

a 

% = 

2.343 k + 

0.4153 - 3.768 

= 

7.15^ k + 

0.1751 - 1.059 Cij^ 

a 6 = 

2.379 k + 

0.05595 + 0.3703 

Th© 

n— 1 or 5th determinant for this 


\ (A25) 


+ a^ag) - (a^a^ - a^ag) ] a^ 

ag(a^a2-aoa3)(a3^-aj^ag) + aiaj^g(a^a^-aoag)-a/a 


3 „ 2 

'6 


(A26) 


The stability boundary can be determined by plotting this determinant 
equated to zero and the coefficient ag of constant or zero-order term 
equated to zero. For this particular example the detailed calculations 
cannot be conveniently 'presented because of their length but stability 
boundaries are presented in figure 1. 

It is customary to restrict the range of values of the variables to 
a practical or sensible limit as has been done in figure 1. This often 
results in part of the bounded stability region being omitted frm the 
chart. Such, an omission is detrimental only from an academic point of 
view. 


After plotting the loci, it is not always Immediately evident which 
bounded region is the stable one. Perhaps the most direct way to 
resolve this problem is to select a point, substitute the values of the 
coordinates into the determinant and the expression for ag and check 
to see if they are positive. Since the point corresponding to the 
values Cju^ = — O.639 and k = 1 is known to lie in the stable region 
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by virtue of the previous calculation, it is apparent that in this case 
the stable region lies between the locus of the upper branch of the 5th 
determinant and the locus of the coefficient ag. The latter is the 
boundary of the stable aperiodic roots and the former is the boundary 
of the stable oscillatory roots. For the value of C of -O.639 used 

in the application of the Hurwitz determinants and the calculation of the 
roots, it is seen that stability exists for control gearings between the 
critical values of O.O7 to 2.3 approximately. These critical values of 

k shall be designated as (k_™) and (k«y,) , respectively. 

low high 


The values of cu giving the frequency spectrum along the oscilla- 
tory boundary are not easily obtained if the technique Just described is 
used for determining the boundaries. In fact, additional calculations 
of the roots must-be made for points lying along the boundary in order 
that the frequencies may be detennined from the imaginary parts of the 
roots. This is a tedious process for 4 th or higher order systems. 

The alternative method for calculating the oscillatory stability 

boundary , also described in the text of this report, provides the fre 

quency spectrum as well as the boundary. In this application substitu- 
tion of the pure imaginary root lo) into the characteristic equation 
yields the following real and imaginary parts which are equated to zero; 

+ agco"^ - a^o)^ + ® 

a^^o)^ — aspfi + ag = 0 

Since Bq, a^ . ag are linear in C„ and k, equations (A27) are' 
inear in and k. Thus it is easy to solve these two equations 

simultaneously for values of Cn^ and k as the numerical value of the 
parameter oj is varied. In figure 1 , it is seen that the range of 
values of cd lying on the part of the oscillatory boundary shown is from 
0 to approximately 22. 



It should be pointed out that the orientation of the static sta- 
bility boundary shown in figure 1 indicates that for a given control 
gearing a decrease in the static stability of the aircraft— autopilot 
combination takes place as -Cjn^, the static stability parameter of the 

aircraft itself, is increased. This undesirable effect plus the fact 
that the airplane by itself had a slightly divergent aperiodic- mode (for 
^1% = —0*639) Indicates that the airplane, as represented by the aero- 
d;pmmic data used, was flying in a critical Mach number range. Had the 
airplane been flying in a subcritical Mach number range, the static 
stability boundary for the aircraft— autopilot combination would have 
indicated Increasing stability as was Increased. 

It is further rioted that part of the stability boundary lies in the 
region where Cjq^ is positive, that is, where the airplane itself is 
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statically unstable. This indicates that an automatic pilot can be used 
to extend the inherent stability limits of the aircraft. 

This method of obtaining the oscillatory stability boundary is more 
expedient the use of the Routh or Hurwitz criterion, at least for 

systems of fourth or higher order. It would seem, therefore, that the 
application of this method plus the use of the locus, for the coefficient 
of the zero-order term is the most convenient way of obtaining this type 
of a stability boundary chart. 


Transient-Response. Analysis 


Although it is possible to obtain an analytical solution for the 
transient response (time history) of the motion of this system following 
the application of a step input in 9 by the operator methods used in 
references ^7 and U8, this solution can be obtained more quickly and 
easily from an anologue computer. In figure 2, the time history as 
obtained from a Reeves Electronic Analogue Computer is presented for the 
case where k=l and Cn^ = —0.639* 

From figure 2, it can be seen that although there is no initial 
overshoot there is an oscillation in the transient before the input sig— ' 
nal is reached, a condition which commonly occurs in higher^ordered 
systems. The steady— state error s^^ is approximately O.O6 or 
6 percent of the step input magnitude. Approximate values for the long- 
period damping constant and frequency were measured and found to be 

( = 2.00 

L 

0J7 = = 10.5 radlans/sec 

‘ 0.60 


which yields 


=10.7 radians/sec 

i j = 0.19 

ti/2 = °*35 sec 


Also by considering the time to reach the first peak in the response to 
be the response time, the speed of response is approximately 2.85 per 
second for a step input of 1°. In a linear system the speed of response 
is proportional to the magnitude of the disturbance. Thus the relative 
rates of response between systems should be measured on the basis of 
units per second per xmit disturbance. 
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Resonance Plot 


In order to plot the frequency-response characteristics of the 
closed— loop hehavior of the aircraft-autopilot combination, it is con- 
venient to put the transfer functions of the aircraft and autopilot into 
forms that can be readily combined to give the frequency-response func- 
tion developed in the text, namely, , - 

f- (io^) =1^ (A28) 

t#i 

Starting with equations (AU), (A5), and (AlO) the operator form of the 
transfer function for the aircraft expressed in determinants becomes 


e 


(-5.325 - 1 . 93 i^D) 

(0.155) 

(-0.639 - O.OI65D) 


0 (5.39) 

0 (0.553 + 1.93^^ D.) 

(0.6299) (0.153) 


(-5.325 - 1.9340) (1.9340) (5-390) 

(0.155) (0.069) (0.553+1.9340) 

(-0.639 - 0.01650), (-0.0340 - 0.01150^) (0.153) 


(A29) 


which reduces to 

i = Ag(0) = ~ -(2.3560^ + 7.1620 + 2.^8^) 

® , 0.0430^ + 0.3180® + 2.8390^ + 0.8730 - 0.178 


(A30) 


Substituting io) for 0 to obtain the frequency response form of the 
transfer function yields . . 


i = Ae(i(o) = - [(-g.356 + 2.383) i (7.162 0))] ' 

^ (0.043 -2.839 <d^-0.181)+1 (-0.318033 + 0.873 03 ) 

+ ibg 
' + iCg 


(A31) 


This function can be written in exponential form. as 


I = A0(ia3) = R^e^^e 


(A32) 
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where 



(A33) 



(hgc 1 — hxCg) 
(hlCx + 1>2C2) 


(A34) 


A corresponding derivation for the autopilot transfer function 
“begins with eq.uatlon (A12) using k = 1. 

(0.058 + D + 37.5) 8 = -37.5 E (A35) 


The operator form of the transfer function is 


| = Ap(D) ^ 


(0.058 d 2 + D + 37.5) 


Substituting ico for D gives 


I - Ap(lu) 


37.5 


(0.058 aP - 37.5 - ico) (di + idg) 


and the exponential form can he written as 


^ - Ap(ico) - Rp6 


where 


^P-= 




e = tan ^ 

^ (di) 


(A36) 

(A37) 

(A38) 

(A39) 

(a 4 o) 


The open— loop transfer function of the combined airplane and autopilot 
can now be written as 




(a 4 i) 
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By virtue of the relationship 


_9 q/E 
9 i • l+(8/E) 


{Ak2 ) 


the frequency-response function of the closed loop becomes 




A( fu) 


Ap(lg)) ^ 

1+A0(io)) Ap(i(o) l+RgRpe^^^®"^^?^ 


(A43) 


or, in the form (A28) desired at the beginning of the development 


where 


and 


9 _ ^2 o^^c 

,9l 


%% 

’/IhE^^^ ^ 2 E 9 EP cos ( Sg+Sp) 


^c = ^9+ ^p 


EgBp sin ( £g+ €p) 
1 +RgRp cos ( £g+ €p) 


(A28) 


(A44) 


(A45) 


Table III shows the calculations carried out in detail from u.' = 0 
to w, =. 20. The plot of the frequency-response characteristics is shown 
in figure 3 * ' 

There is good correspondence between the resonant frequency, 
approximately 11 radians' per second, and the long-period natural' fre- 
quency, 10.64 radians per second, as calculated from the roots of the 
characteristic equation and measured' from the transient-response data. 
The steady-state error is computed (table III) to be 8.2 percent com- 
pared to the less accurate value of 6.0 percent measured on the tran- 
sient response. The relatively high peak amplitude of approximately I .7 
indicates that the damping of the dominant oscillatory mode is not very 
great. This is borne out by the damping ratio of approximately 0,2 for 
the long— period oscillation obtained from the roots and by the transient 
response. The assumed relationship between the peak amplitude and the 
peak overshoot was not borne out in view of the fact that no initial 
overshoot of the transient response occurred, although a relatively higli 
peak amplitude is shown in figure 3 ] The oscillation of the transient 
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before reaching 9 ^, which commonly occurs for systems of third or higher 
order having high peak ajnplitudes, did occvir, however.. The relatively 
high value of the resonant frequency indicates a high speed of response. 


Greenberg's Method 


This method will be applied to calculate the critical control gear- 
ings (kf-), . , and (kcr), 0™ fixed at -O.639. In 

high low ^ 

Greenberg's notation the frequency-response form of the transfer function 
for the aircraft is 


i = A. e-l«r 

8 Kr 

as contrasted with equation (A32) where 


i = 


Also the autopilot response is designated as 

^ . Kpe^'P 

e P 

which corresponds to equation (A38) where 

^ -D 

P 


(a46) 


(Akj) 


The two ^phase angles are not equal, • however . In Greenberg* s method 

the input 0 Is assimed to be zero and the error is assumed to be 
equal to the output, that is, 

E = a (a48) 


rather than the usual assumption that 


or for = 0 


E = -a 


(All) 


E = -9 


(A49) 
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Magranmiatically this arrangement appears as follows: 



and in order for the aircraft-autopilot comhlnation to he stable It Is 
necessary to use the autopilot eq,mtion in the form of (a8) 

( 0 . 058 d^ + D + 37.5) 8 = 37.5 k E 

With the conventional notation and relationships the point — l40i is 
the critical stability point as established by Wyquist. In other words, 
when the phase angles add up to l80°, ® 2 = n the product of the 

magnitudes Eg, Ep, and k must be less than or equal to 1. Thus 
for comparison, and writing Greenberg’s autopilot phase angle as Cp , 

% = ^ (A50) 


^0 = - ^r (A51) 

and 

t ' , • ■ 

Epe^^P =-Kpe^^Pi . (A52) 

or 


so 


Rpe^S = Kpe~^" e^^Pi 


(A53) 


Ep = Kp , (A54) 





(A55) 


— rt 
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In Greenberg's notation, at the point 


^0+ = 


- "r^ "P, - « = 


(A56) 


or 



(A57) 


the magnitudes must satisfy 


k 


Kr. 




(A58) 


for stability. 


To find the critical control gearings, therefore, the resonance 
plots of the autopilot and aircraft are compared and at the frequencies 
where the angles and are equal, the amplitudes Kp and Kj. 

can be determined and substituted into the following expression; 


^r ~ 



(A59) 


Thus if the control gearing is set equal to its critical value, equation 
(A58) is satisfied. 


k ^ ^ 5 p = 1 

cr Kp Kr 

The values of the critical control gearings or gain margins based on the 
values of Kp and Kj. at the critical frequencies shown in figure k 
are 



2.80 

1.26 


2.22 and (kcr). = = O.O75 

low 1.00 


The phase margin, also available on this type of plot, is 50°. The com- 
putations are available in table III. 


Open— Loop Transfer Locus Plot (Wyquist Diagram) 


The transfer locus showing the open— loop response of the aircraft- 
autopilot combination for Cn^ = — O.639 and k=l was obtained by 
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plotting the individual transfer functions (A32) and (A38) and combining 
them vectorially. The results are shown in figures 5 and 6 at different 
scales and the computations are included in table III. Figure 5 shows 
a small-scale plot giving fairly good detail of the transfer locus in the 
region of high frequencies. The phase margin and the higher value of 
gain margin are evident as 48° and 2.22, respectively. From figure 6, 
which shows the transfer locus plotted for o)=0 as well as the high fre- 
quencies, it can be seen that the lower value of critical control gearing 
is 0.076. It should be pointed out that although the airplane auto- 
pilot are combined in the form of a single— loop system, the airplane by 
itself is unstable and should be considered as a multiple— loop system 
within the single loop. Thus the check of stability on the transfer 
locus plot should be made using the generalized Nyquist criterion relat- 
ing the encirclements of the — 1+Oi point to the number of zeros minus 
the number of poles as described in reference 11. 

The effect of the simplifications in the equations of motion, used 
in the section of the report discussing damping, upon the frequency- 
response characteristics was of sufficient interest to warrant a compar- 
ison. The simplified transfer function of the aircraft in operator form 
is 


Aq (d) = 0.6299 (3.323 + 1.933 D) 

® ' (o.967)(0.0231) d (d^ + 7.12 D + 63.4) 

The frequency response form Is 

Aq (ico) + 1.218 

^ — 0.159 + (— 0.0223 + 1 . 4 l 6 ) icD 

_ e 1 + i e2 
fi + if 2 

In polar form. 


(A60) 


(A6l) 


(A62) 


A0^(1o3) = Kq^ 


(eifi+esfa) +(®2f 1 — eif2)^ 


+ fc 


ei tan~i 


(egfi-eifg) 

(exfi+e2f2) 


A0g(iai) = 


(A63) 
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The calculations ,£ire carried out in table IV. The simplified airplane 
transfer locus is plotted eind combined with the , autopilot transfer . locus 
in figure 7* Comparison of these results with those shown in figure 5 
indicates that the simplification was noticeable only at the low fre- 
quencies. No check on the lower critical control gearing was possible 
but the values of the phase margin and the higher critical control 
gearings are virtually the same. 


Logarithmic Plots 


The commonly used logarithmic form for plotting the transfer locus 
has the decibel as the unit of magnitude. To get the magnitude of the 
transfer function i«ito this form the logarithm to the base 10 is taken 




CO 


and multiplied by 20. Thus the log modulus is defined, in decibels, as 

Lm I A(icn) | = 20 logio lA(icu) | (A64) 


> The phase angles are not altered. 

For the aircraft transfer function 


I = AgCico) = 


Id + ilD2 


Cl + iC2 

■bhe log modulus is 

Lm I Ag(ico) I = 20 log^Q -v/bi^+bs^ - 20 logio v^Ci^ + 

= 10 logj^Q (bi^+b2^)-10 logj^Q (ci2+ t^) 


(A 3 I) 


\ (A 65 ) 


Writing the phase angle in terms of the numerator and denominator com- 
ponents yields 


— 1 bg 


Ang (bi + ibg) = tan — 

bi 

Ang (ci + icg) ^ = tan”^ ^ ^ 


> 


(a66) 
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For the autopilot transfer fimction (A37) . 


I ,= Ap(to) 

E ^ di + ids 


the log modulus is 

Lm I Ap(ia>) I = 20 log^o dg - 20 log^^g 4 di^ + dg^ 1 


= 10 logj^o 


V 2 . ^ 2 / 


and the phase angle is 


Ang { ') = tan ^ f ') 

^ \ di + ids/ \ di / 

The over-all open-loop transfer fimction is 

9 - ' 

- = A0p(i(o) = Aq (ioo) Ap(io») 

The log modulus, therefore, is 

Lm I Agp(i(o) I = Lm | Ag(io)) | + Lm | Ap(io)) | 

= 10 lOgio (tl^+t2^)+10 lOgjo (C 1^402^) + [ 


and the phase angle is 


Ang A^p(ia>) = tan ^ ^ + tan ^ ^ ^ + tan ^ ^ 


(A67) 


(a67) 


(a68) 


(A69) 


(A70) 


(A71) 
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The calculations are given in table V. The log modulus and phase- 
angle ciorves plotted , against the log of the frequency are shown in 
figvire 8. A graph showing the log modulus versus the phase angle is 
also presented as figure 9* From this curve the gain margins and phase 
margin are seen to be 6.9 decibels below the zero level, 22.4 decibels 
above the zero level, and 48°, respectively. The value of the gain 
margins, 6.9 £^iid —22.4 decibels, correspond to the values of 
*^™hi h “ 2.22 and Ompow ” 0*076 which were found from the polar plots. 


Approximate Factorization of Quartic 


As mentioned in the body of this report, the log modulus analysis 
is particularly convenient if the system being analyzed is composed of 
only first— and second— order transfer elements. To take advantage of 
this situation in the case of longitudinal motion, the denominator of 
the aircraft transfer function often can be approximately factored into 
two quadratics by the method used in reference 39* 

If the quartic is of the form 

+ aiD^ + + a 3 D + = 0 (A72) 


it can be written as 

(D^ + ZjD+za) (D^ + Z 3 D + Z 4 ) = 0 


(A 73 ) 


where 


Eq = 1 

El = Zi + Z3 

8-2 = Z3 -f Z2 + Z4 

E3 = ZiZ^ + Z3 


(A74) 


In reference 43, the following approximations are recommended for all 
practical cases of conventional’ aircraft: 
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zi = ai 

Z2 = as 

rr _ ®'3 

Z3 ^ 

a .2 0-2 



as 


From eqmtion (A30) the denominator of the aircraft transfer func- 
tion is 

0.01^31 (d'* + l.h 2 h D® + 65.839 + 20.235 D - 4.190) (A76) 


Substituting in the approximate relationships yields 

y 0.0431 (D^ + 7.424 p + 65.839) (D^ + 0.314 D - 0.064) (A77) 


Thus in place of Cj + ics in equation (A31)j the following 
quadratic factors can be used: 

(g-1 + igs) (^1 + ihg) = (Zi + ± 1 ^) 

where with ico substituted for D 

g^ = 2.838 - 0.0431 0^ 
gg = 0.3320 Oi 
h^ = -(0.064 + (a) 
hg = 0.314 CD 

The calculations are given in table VI and the curves shown in figure 10. 
A comparison of the plots of the log modulus of and + iCg 

can be made using figures 8 and 10. It is seen that in- this particular 
case the approximation is very good. 



) (A75) 
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APPENDIX B 


SYMBOLS AND COEFFICIENTS 


A0p 


A, Aj^, As, A3 arbitrary transfer functions 

Ac closed— loop transfer function for aircraft and autopilot 

combined 

Ap autopilot transfer function 

Aq aircraft transfer function for longitudinal motion 

open— loop transfer function, the product of Aq and Ap 

simplified aircraft transfer function for longitudinal 
motion 

amplitude ratio (reference 9) 
coefficient In arbitrary transfer ftmction A 
coefficient in arbitrary transfer function A 
drag coefficient 


^0 

AP 

B 

C 

Cd 

gl 

Cm 


8 




a 


lift coefficient 




pitching-moment coefficient 
thrust coefficient 


^ moment y 
\ qSc / 


^ thrust A 

\ qS / 


rate of change of lift coefficient with angle of 
attack 


\ 


.rate of change of drag coefficient with angle of 

)aj 


/acD'' 

attack I 

\ oc 



rate of change of pitching-moment coefficient with 
angle of attack 



rate of change of pitching-moment coefficient with 
pitching velocity 




differential operator 


error signal {9^S) 

arbitrary transfer function (reference 12) 

frequency dependent part of a transfer function 
(reference 3) 

gain margin 

' j 

moment of inertia 


moment of inertia of motor 


frequency invariant part of a transfer fimction 
(reference 3)j also the loop sensitivity or gain 


nondimens ional moment of inertia parameter 


Ky^pS 

me 


aircraft amplitude ratio (reference 7 ) 
autopilot amplitude ratio (reference 7 ) 


radius of gyration of aircraft about axis of pitch 
Laplace transform of [] 
log modulus 

constant magnitude parameter for closed-loop response 
constant phase angle parameter for closed— loop response 


period, seconds 
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PA phase angle (reference 9 ) 

^ performance operator (reference 9) 

E, El, Eg, Es, E4 arbitrary magnitudes 


% 

«^s 

S 

Ss 

St 

T, Ti, T2». 
V 

y 


amplitude of autopilot transfer function in polar form 

amplitude of aircraft transfer function in polar form 

amplitude of simplified aircraft transfer function in 
polar form 

wing area 

torque— speed sensitivity of motor 
torque— voltage sensitivity of motor 
arbitrary time constants 
velocity of aircraft 

arbitrary transfer function (reference 11) 


®-0^. 


®' 4 J ®' 5 ^ 

' coefficients of characteristic equation 

®-n 

' . ' 

bi, ba 

coefficients in Ag. 

c 

mean aerodynamic chord 

Cl, ca 

coefficients in Ag 

d 

logarithmic decrement 

CVI 

H 

coefficients in Ap 

• 

e 

2.71828 

ei, ea 

coefficients in Aa 
03 

f 

damping coefficient (viscous) 

^cr 

critical damping coefficient 
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damping coefficient of motor 


coefficients in 
acceleration due to gravity 

coefficients in an approximate factorization of 
stability quartic for aircraft 

coefficients in an approximate factorization of 
statillty quartlc for aircraft 




proportional control constant 


control gearing 


6(i(o) 

<^0 


critical control gearing 


derivative control constant 


Salii-«2^2 

ggh^+gihg 


mass of aircraft 


arbitrary quantity 

a variable introduced in the laplace transformation 


free— stream dynanlic pressure 




roots of a second— order equation 
steady— state error 
amplitude of initial overshoot 
time^ seconds 


time to damp to half amplitude 
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tr 

u 

u* 

X 


ZSj Z3, Z4 


6 

6a 

6 e 

8r 


response time 

longitudinal perturtation velocity 
nondlmensional longitudinal perturbation velocity 

arbitrary dependent variable 

coefficients in approximate factorization of aircraft 
stability quartic 

angle of attack, radians 

phase margin, degrees 

arbitrary control deflection, radians 
aileron deflection, radians 
elevator deflection, radians 
rudder deflection, radians 
arbitrary phase angle, degrees 



%• 

^08 


5 

0 


phase angle for aircraft-autopilot closed— loop transfer 
function, degrees 

phase angle of autopilot response, degrees 

phase angle of aircraft response (reference 7), degrees 

phase angle of aircraft response, degrees 

phase angle of simplified aircraft response, degrees 

t 

f f 

damping ratio ( 

' ^cr 

output angle of pitch, radians, unless otherwise noted 
input angle of pitch, radians, unless otherwise noted 
arbitrary exponential coefficient 
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P 


T 


'"d 

9 

<p(t) 

(U 


roots of a fifth-order characteristic equation 

spring constant 
density of air 




aerodynamic time 


time lag 


angle of hank 


arbitrary function of time 


seconds 


angle of yaw 

actual frequency ^ 1 —^*" J, radians per second 

natural frequency, radians per second 




Subscripts 


I 

o 


long— period oscillatory mode 
initial ralue 


/ 


s 


short— period oscillatory mode 
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KESONAWCE PLOT CALCULATIONS 
[See equations (A28) through (a 45) and fig. 3.1 
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CALCULATIONS FOE THE SIMPLIFIED TRANSFER LOCUS OF TEE AIRCRAFT 
,[ See equation (a 63) and fig. 7. ] 


NACA TN 2275 



CM 

CM 





^ \0 lf\ CO^ COHCUVOVO^HVOO 
- 4 - D-irNH-itONHirsCrsiHro 0 \C\JI^ 


r-l 

CM 


0 

0 

0 iH 0 \^ 0 H 0 0 l>-CO ITN 0 -=t CVJ 

0 

1 

+ 

jfop-:±OJ<^CVjd^ 

1 1 1 1 1 1 1 1 


H 

rH 



ch 



(M 

H 



0 

0 

' 



CM 




ON OMf\ H -:t CVJ OO IPv O CVJ C— 

00 VO OJ 0 H VO 0 \-:t 0 VO CJ C\J 



H/ 


0 

, OOCOir\-=f- 4 ’VOVOHON c— ^ 00 0 ON 

© 



• o************ 

0 OJ ^ l^- 0 t-vo CO cnco 0 >VO H ITN 
II HHCy/C^lOt-ONCVlVO 

1 1 1 II 

H 

CM 


0 




CM 




t>-VO VO OJ C\J \r\J=t H ^ yo ^ iT\ CT) 



CM 

ON VO -It cn IfN 0 CVJ CO ITNVO ^-VO VO 


0 

CVJ HirNOOtfNHOOVOCOOVO-^OO OJ 






+ 

H 

0 HCVJ-:tOO(^JCVJt^VOirNHCVJVOCJN 
H H T cpop ITN 0 VO 

1 1 rH H CO 



H CO 



H 

1 1 


0 



CM 

C 7 N ON VO CO 00 VO 00 OJ 0 0 ON ON 0 



CM 

ON VO OJ-it OOH-cj* OnHVOOOVO OO 



ON 0 O-=J- OO-t OJ C\i CO LTN H’ 0 O- 

© 

+ 

0 H -It VO ITN COCO ON 0000 



H CM OO-it VO 0 H ON VO 

CM 

iH rH 00 iH LTN 



H 

•V *\ 



-It VO 




CVJ 




lTN- 4 - ON-:t IfNVO ITNCO H VO O O ^ O 
0 C 7 N-it lT\-:t 00 ON C^VO CTn H 
t^oOOVOVOCViOJVOCVJOlTNHOO 

© 

CM 

• ••••••••••-• ** 

0 H OJ OJ oo-it -It 00 OJ 1 cncp 0 

'1 




CJN ON t^vo H -f ITN-It H VO ( 0 \ 0 0 0 
00 ITN lA 00 OO-it t^CVJ ONl^t^OOO 

0 

H * 

0 H mVO-^ iTNONO-t^-HCO ONl^VO 

Ol 1 1 

1 1 II I ^ 







ONCO l^-VO -t- CVJ 0 00 VO CVJ 0 0 0 
OHOJ OOITNO-ONO CVJ-^VOOO l^'VO 

©) 



VOCVJOO-t-VOOOOOOlTNC^ONHCVJ CO 

0 

• m % m m m m 

Ky 



0 iH H CVJ 00-it VO C-“CO ON 0 OJ 00 -It 

* 



iH rH iH CVJ 

© 



ITN -ITN 

2 


0 H H CVJ CO-t ITNVO 1^00 CJnO lA 0 

»H iH OJ 



mck TN 2275 




CO 

CD 


00 000000 000000 

VO CO H CT\CO ONLPv COCA H CVJ ^-CO 

VO ovoiAir\c-irMT\ d-oo vo o -=i- h 


^OOooOOOHONHVOOJ D-OCOON 
Q? i^r* ^rVQ CY-t^CAOOJ couaooh 

I I IT 


CO 

CD 


CVIHOOJ CAOmO 
VO-=l--4 -OCOOJC\JOOO CA-:t O l>-(^ 
00 LT\ CO GO CU CVJ OJ CVi CvJ O CAOO CO CVJ 

^ tPsCVJHHHHHHHH 




00 000000 000000 
0-=fHO-:t-:t-=J--=J-Oir\VOOVOON 
ON lA CO ir\ Ov ITS VO ON CO LTN CVJ m ooco 

* ••••••••• •o«* 

CVI VOH t>~COir\HH LTvCACVJCVJVO-^J- 
VO LA Lf\-=J- -4- lAVO ^-co O H Co3- 
Cjj C|JC|JCjJCjJCVJCjJCjJCjJCVJ(p 


{o\ 

I! 

CQ 

CD 

w 


A^HOHVOO-!>-H 
'^'^d:OI>~CACOCOCOH AOCOA 
CO AGO COCVJHHHHOOOC--CVJH 

• *****o***o««« 

? A C\J H H H H H H H H 




A O 

O HAOCOHOOO 

OOO^-itCVJt^D-^vO ^-A 

OCJ AD-cOACAO-^-^^OJlA 


HHCVJ-^AO-O ACVJ 


H H 


O 

O 


CVJ cvi o 

•\ 

H 


CM 

© 

+ 

(VI 

&] 


VO A i>- l>- 

AVOOOOCVJACVJVOA 

!>-OHA00H-=J*OO 

O A ACOtr-CO-:^ ACOCVI OVO CO 
CUVOCVJ O^AAACVJCO HCVj 
H COCO VO A A A-:t A 


CVJ A O 
H 


CO O 
CVJ A 


o o 

t-o 

CO O 

•V *N 

O A 
^ VO 

^ o 

•\ *N 

H VO 

iH 


A A 

H H CVJ. CO-^ AVO D~CO 0\ O A O 

r—J rH CVJ 














KACA TN 2275 


97 


TABLE V 


CAKIULATIONS FOR THE LOGARITHMIC PLOTS 
[ See equations (A 64 ) through (A71) and fig. 8. ] 



© 

© 

® 


® 


bi^+hs^ 

C12+C22 

di^+d. 2 ^ 


10 logio© 

di^+d^^ 

0 

5.67464 

0.03263 

1,406.250 

1 

0 

.5 

16.03440 

.94465 

1 , 405.413 

1.00059 

.0026 

1 

51.27673 

9.16085 

1,402.903 

1.00238 

.0103 

1.5 

123.89144 

40.35783 

1,398.730 

1.00538 

.0233 

2 

254.69419 

118.26430 

1,392.904 

1.00958 

.0415 

3 

815.76008 

530.65061 

1,376.372 

1.02170 

.0932 

4 

2,067.52987 

1,483.1020 

1,353.511 

1.03896 

.1660 

5 

4,476.28100 

3 , 224.141 

1,324.602 

1.06163 

.2596 

6 

8 , 641.515 

6 , 245.461 

1,290.010 

1.09010 

.3747 . 

7 

15,295.527 

12,028 

1,250.177 

1.12484 

.5108 

8 

25,305.667 

24,647 

1,205.629 

1.16640 

.6685 

9 

39,669.360 

53,031 

1,156.971 

1.21545 

.8472 

10 

59,519.605 

118,582 

1,104.890 

1.27275 

1.0473 

15 

290,032 

3,545,111 

822.802 

1.70909 

2.3277 

20 

904,174 

39,669,789 

604.490 

2.32634 

3.6667 
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TABLE V.- COMTHDED 





© 



(A) 

10 logio 0 

0 

0 

0 

© 


tan 1 

— C2/ Cl 

0 

7.U3 , 

14.864 

0 

-180° 

0 

.5 

12.049 

0.248 

1.99669 

-116. 60° 

.44636 

1 

17.098 

-9.619 

275.29600 

-90.21° 

.18557 

1.5 

20.968 

-16.058 

-3.67979 

-74.79° 

.03594 

2 

24.062 

-20.726 

-e. 03371 

-63.81° 

-.07523 

3 . 

29.155 

- 27.248 

- 1.14132 

- 48 . 10 ° 

-.27099 

k 

33.155 

-31. 711 

-.81108 

-39.05° 

-. 49182 

5 

36.508 

- 35.084 

■ -.63348 

-32.35° 

-.80662 

6 

39.365 

- 37.948 

-.52118 

-27.53° 

-1.37452 

7 

41.873 

-40.803 

-.44333 

-23.91° 

-2.90004 

8 

44.032 

-43.918 

-. 38601 

- 21 . 11 ° 

29.85025 

9 

45.984 ■ 

-47.294 

-.34196 

-18.88° 

4.27415 

10 

47 . 746 

-50. 741 

-.30703 

-17.06° 

2.11691 

15 

54.623 

— 65.496 

-.20353 

- 11 . 44 ° 

.68800 

20 

59. 731 ' 

- 75.984 

-.15235 ■ 

- 8 . 6 ^ 

.44138 
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TABLE CONCLDliED 


1 ■ 

0 


@ 

0 

0 

0) 


-is/dl 

- tan“^ @ 



0 

-180° 

0 

-180° 

22.307 

0 

0 

.5 

-155.95° 

-.01333 

-180.76° 

12.299 

-93.31° 

1 

-169.^9° 

-.02670 

-181.54° 

7.489 

- 81 . 24 ° 

1.5 

-177.93° 

-.04013 

-182.30° 

4.933 

- 75 . 02 ° 

2 

- 184 . 30 ° 

-.05366 

-183.09° 

3.378 

-71.20° 

3 ' 

-195.15° 

-.08112 

- 184 . 64 ° 

2.000 

-67.89° 

k 

-206.20° 

-.10937 

-186.25° 

1.610 

. -71.50° 

5 

- 518 . 88 ° 

-.13869 

-187.91° 

1.684 

- 79 . 14 ° 

6 

r 

-233.96° 

-.16943 

-189.63° 

1.792 , 

-91.12° 

7 

-250.96° 

-.20197 

-191.43° 

1.581 

-106.30° 

8 

-268.90° 

-.23677 

-193.34° 

.783 

-123.35° 

9 

-283.16° 

-.27437 

-195.35° 

-.463 

-137.39° 

10 

-295.29° 

-.31545 

-197.52° 

- 1.948 

-149.87° 

15 

-325.47° 

-.61349 

-211.51° 

-8.545 

- 188 . 42 ° 

20 

-336.19° 

-1.39860 

-234.29° 

-12.586 

-219.14° 
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© 

© 

© 

© 

0) 

gx 

S2 

111 

0 

2.838 

0 

-0.064 

.5 

2.827 

.166 

-.314 

1 

2.795 

.332 

-1.064 

1.5 

2.741 

.498 

-2.314 

2 

2.666 

.664 

-4.064 

3 

2.450 

.996 

-9.064 

4 

2.148 

1.328 

-l6.o64 

5 

1.761 

1.660 

-25.064 

6 

1.286 

1.992 

-36.064 

7 

.726 

2.324 

-49.064 

8 

.079 

2.656 

-64.064 

9 

-.653 

2.988 

-81.064 

10 

-1.472 

3.320 

-100.064 

15 

-6.859 

4.980 

-225.064 

20 

-14.402 

‘ 6.640 

-400.064 


© © ® 


ga/gi lieAi 


_i —Es 


0 

.157 

.31^ 

.471 

.628 

.942 

1.256 

1.570 

1.884 

2.198 

2.512 

2.826 

3.i4o 

4.710 

6.280 


0 

.059 

.119 

.182 

.249 

.407 

.618 

.943 

1.549 

3.049 

33.620 

-4.576 

-2.255 

-.726 

-.461 


0 

-.500 

-.295 

-.204 

-.154 

-.104 

-.078 

-.063 

-.052 

-.045 

-v039 

-.035 

-.031 

-.021 

-.016 


-13.79° 

-22.15° 

-31.72° 

-43.34° 

-57.14° 

-71.85° 

-88.30° 

-102.33° 

-113.20° 

-144.03° 

-155.25° 


-180° 

-153.45° 

-163.56° 

-168.47° 

-171.26° 

-174.06° 

-175.54° 

-176.38° 

-177.03° 

-177.42° 

-177.76° 

-178.00° 

-178.21° 

-178.80° 

-179.09^ 
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Static stability factor f-Cm^ 


Figure A— Stability boundary chart for the aircraft-autopilot 
combination. 
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f^crJ,ow= 0.076 


\mm 


Frequency^ (if 

Figure 4.- The critical control gearing ^ as (Stained from the 
frequency-response characteristics of the aircraft and autopilot. 
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Figure S.^The open-loop transfer tool of the aircraft , the 
autopilot, and the aircrafhautopiiot combination. 



Figures.- The open-loop transfer locus of the aircraft-autopilot 
combination showing the low frequency region. 
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Figure 7.- The open-loop transfer locus of the aircraft using 
simplified response equations, and the locus of the 
corresponding aircraft-autopilot combination. 
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Frequency, oj 



Figure 8.- Logarithmic plots of the transfer loci of the autopilot d^/(d,4idg), 
the aircraft (bi+lt^l/fci-rlc^), and the combination representing the 
controlled airplane. 
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showing the gain and phase margins. 
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Figure 10.- Logarithmic plots of the transfer loci representing the 
approximate factorization of the aircraft stability quartic 
c,-hiCg into the quadratic factors (gi+igg)(hi+ih^)sl^+ii^. 
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